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Abstract 

Starting from the famous Pauli problem on the possibiUty to associate 
quantum states with probabihties, the formulation of quantum mechanics 
in which quantum states are described by fair probability distributions 
(tomograms, i.e. tomographic probabilities) is reviewed in a pedagogi- 
cal style. The relation between the quantum state description and the 
classical state description is elucidated. The difference of those sets of 
tomograms is described by inequalities equivalent to a complete set of 
uncertainty relations for the quantum domain and to nonnegativity of 
probability density on phase space in the classical domain. Intersection 
of such sets is studied. The mathematical mechanism which allows to 
construct different kinds of tomographic probabilities like symplectic to- 
mograms, spin tomograms, photon number tomograms, etc., is clarified 
and a connection with abstract Hilbert space properties is established. 
Superposition rule and uncertainty relations in terms of probabilities as 
well as quantum basic equation like quantum evolution and energy spectra 
equations are given in explicit form. A method to check experimentally 
uncertainty relations is suggested using optical tomograms. Entanglement 
phenomena and the connection with semigroups acting on simplexes are 
studied in detail for spin states in the case of two qubits. The star-product 
formalism is associated with the tomographic probability formulation of 
quantum mechanics. 
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1 Introduction 



Pure quantum states are usually associated with wave functions [T] or vectors in 
a Hilbert space [2] . Mixed quantum states are associated with density matrices 
[3] or density states [3]. Pauli |5ii6J posed the problem whether it was possible to 
associate quantum states with probability distributions as it happens in classical 
statistical mechanics. Pauli problem was more concrete, namely, is it possible 
to reconstruct the quantum state (i.e., the wave function) from the knowledge 
of the probability distribution for the position and the probability distribution 
for the momentum? The answer to this particular question is negative (see 
the discussion in Reichenbach's book J7, a-nd a recent example in [8J). But the 
general idea of Pauli to associate quantum states with probability distributions 
was implemented by introducing the tomographic probability representation of 
quantum states [9j. This representation is based on the Radon transform [lOj 
of Wigner function [TT] , suggested in [T^l [13] to connect the measurable optical 
tomographic probability [TU [T5] to reconstruct the Wigner function of a photon 
quantum state. The mathematical nature of the tomographic probability rep- 
resentation was clarified in [HI [H] (but see also, e.g., [TO l ) . 
In quantum mechanics we have the conventional Heisenberg and Schrodinger 
representations. The tomographic probability representation is another one. 
The physical properties of quantum systems can be studied in the tomographic 
probability representation as well as in the Heisenberg and Schrodinger represen- 
tations or in the Feynman representation [24, based on the use of path integral 
as main ingredient of the quantum picture. Since the tomographic picture deals 
with probabilities which describe the quantum states, we will show how impor- 
tant quantum aspects as uncertainty relations and superposition principle can 
be described in terms of tomographic probabilities. We suggest a method of 
checking the Heisenberg uncertainty relations using quantum state tomograms. 
It is worthy to mention that all the available representations of quantum me- 
chanics are equivalent (see, e.g., the review [25_). One cannot say that some 
representation is better or worse than the others. Nevertheless, each repre- 
sentation has peculiar properties due to which some quantum aspects become 
clearer and simpler than in other representations. For sure, the superposition 
principle can be formulated in the easiest and clearest form using the natural 
linear structure of a Hilbert space whose vectors are realized by complex wave 
functions. The tomographic probability picture is very natural for problems of 
quantum information and quantum entanglement, so each picture has its own 
merits. 

We shall try to present in a pedagogical style the construction of the to- 
mographic probability representation both for discrete (spin, qubit) variables, 
studied in [26[ 127] , and for continuous variables like position and momentum 
following [HI [m [171 [2H1 [29l [30] . The tomographic probability can be used also 
in classical statistical mechanics [3TJ [32]. The tomogram of a classical state 
is the Radon transform of the standard probability density on the classical 
phase space. In this setting both classical and quantum states can be described 
by tomographic probability distributions. The difference between classical and 
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quantum states in the tomographic description is related to the diflFerent phys- 
ical constrains the state tomograms have to satisfy in order to be either in 
the classical or in the quantum domain. In the tomographic representation 
the quantum-classical relation is formulated in terms of properties of the tomo- 
graphic probability densities. It will be obvious that the ambient space of Radon 
transforms of functions defined on the phase space contains the subset of quan- 
tum tomograms (tomograms admissible in the quantum domain), the subset of 
classical tomograms (those admissible in the classical domain) and the subset of 
those which are not admissible neither in quantum nor in classical domain. The 
subsets of quantum and classical tomograms have a not empty intersection. In 
this work we study the properties of tomograms providing a characterization of 
classical and quantum domains as well as of their intersection. 

The paper is organized as follows. After a preliminary section 2 in which 
the Pauli problem is considered and its tomographic solution is introduced, 
there are two main parts. Part one is devoted to discuss quantum mechan- 
ics formulations on phase space and their relations to tomography. It contains 
nine sections. Weyl systems are discussed in section 3.1. Wigner functions are 
considered in section 3.2, and their transformation properties under the action 
of the automorphisms of the Weyl-Heisenberg group are discussed in section 
3.3. Tomograms and Radon transform are studied in section 3.4. Classical and 
quantum probability distributions are considered in section 3.5. State recon- 
struction procedure is studied in section 3.6. Tomographic families of rank-one 
projectors are introduced in section 3.7, while a general and abstract setting 
of tomographic maps is presented in section 3.8. Finally, a unified approach 
to construct the most commonly used tomographic families of observables is 
given in section 3.9. In the second part we describe quantum mechanics in the 
tomographic picture. It contains six sections. In the sections 4.1 and 4.2, the 
superposition rule and the uncertainty relations respectively are discussed in 
the tomographic probability representation. In the section 4.3 some examples 
of classical or quantum distributions are presented. In section 4.4 we study 
basic equations (time evolution and eigenvalue equations) in tomographic rep- 
resentation. An application to entanglement and separability on examples of 
two qubit states is given in section 4.5. Bell inequalities and association with a 
semigroup structure are considered in section 4.6. In section 5 some conclusions 
and perspectives are finally drawn. 

2 The Pauli problem: the original formulation 
and the tomographic solution 

In his book [6], after having introduced the wave function in the position rep- 
resentation by means of tjj {x) and in momentum representation by means of 
(p {pj , along with the probability densities 

W{x)^\ij{x)\' ^^PixYijix); Wip) = ^{p)*^{p), (1) 
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Pauli claims: These functions ip (x) and fip), which are usually called 'prob- 
ability amplitudes' are not, however, directly observable with regard to their 
phases; this holds only for the probability densities W (x) and W {p). In a 
footnote Pauli states: The mathematical problem as to whether for given func- 
tions W (x) and W {p), the wave function V', if such a function exists, is always 
uniquely determined (i.e. if W (x) and W [p) are physically compatible) has 
still not being investigated in all its generality. [6] 

As it was shown in Reichenbach's book [7] , the negative answer to the origi- 
nal Pauli problem may be given by means of counterexamples. For instance [5], 
consider the two squeezed states with wave functions in the position represen- 
tation {h = 1): 



ijjiiq) = NeyiY>{-aq^ +i(3q), 



(2) 



V'2(g) 



Rea > 



7Vexp(-Q!*g^ -I- ifiq), 



f3 = f3* ■ N 



a + a* 



TT 



and in the momentum representation 



exp 



N 



: exp 



4a 

_{J3±pf 
4a* 



(3) 



One then has 



\M<i)\' 



|7^2(g)|' = |7V|%xp[-(a-Ha*)g2] 



(4) 



\Mpt = \Mp)\' = 



\N\ 



a 



■ exp 



-W + P) 



Aaa* 



The fidelity I (Vi|V;2)r 



/ = 



2^; 



aa' 



One can see that / 7^ 1, which means that the states are different. 

To understand the reasons why the knowledge of the two marginal distri- 
butions of position and momentum are not sufficient for reconstructing a state, 
consider the family of dimensionless observables, depending on two real param- 
eters /i, : 

u)^nQ + vP , (5) 

where, restoring the Planck constant h, Q and P generate the Weyl-Heisenberg 
algebra [(5,P] ~ ihl. The spectrum of X{ii,i^) is the real line, which we 
parametrize by X, with corresponding improper eigenvector \Xiii>). In the 
position representation. 



X{ii, v) — —ihv— 
dq 



(6) 
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and its (improper) eigenfunctions may be chosen as 



y/2nh\i^\ 



(7) 



The normahzation is such that {X'\X) — S{X — X'). Analogously, in the mo- 
mentum representation, the eigenfunctions may be chosen as 



1 



(8) 



Now, define the symplectic tomogram of a (normalized) pure state 1-0) with 
respect to the family X{^,v) as (see, e.g., [55] ) 



T^{x,^Ji,u) = \{x^ly\^)\^^ 



1 



1 



'4'{p)e ^2h^P +^n^Pdp 



2Trh\fi\ 

where the Fourier transform of the wave function has been introduced 

1 



■0(y)e * R dy. 



(9) 



(10) 



In other words, %p{X,fi, v)dX is the marginal probability such that a measure 
in the given state \'ip) of the observable X(/i, z^), with fixed /x, j/, has values in 
{X, X + dX). Of course, as IV') is normahzed. 



j T^{X,ii,v)dX = l 



(11) 



independently of the chosen values of /i, v. 

Then, it is apparent that the Pauli problem amounts to reconstruct the 
given state from the knowledge of the two marginal probability distributions 
of position, 7^(X, 1,0), and momentum, 7^(X, 0, 1). As a matter of fact, the 
following reconstruction formula holds (see, e.g., |34j ) for the density matrix 
p{q,q')=^{q)r{q')-- 

p{q,q') = ^y'r^(X,^,i^)e4^-^(^+'?')]dXdA. (12) 

= / %p{X, ^, i/)e^P(«~'?') exp [i{X - pq- vp)] dXd^idvdp 

(27r)^ J 

Thus, the answer to the Pauli problem is negative because the reconstruction 
requires the knowledge of many different marginal probability distributions, 
corresponding to many different observables of the family X{^,v). A minimal 
set of such observables is called a quorum 35 , and the characterization of a 
tomographic quorum will be discussed in the next sections. 
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3 Quantum mechanics on phase space and to- 
mography 



3.1 Weyl systems 

The notion of Weyl system is useful as a tool to formulate a quantum version 
of classical hamiltonian mechanics by using the symplectic form on a linear 
classical phase space. We will use it in the following as a suitable setting to 
discuss the quantum and classical tomographic maps and their relations. Here 
we briefly recall the definition of Weyl system, mainly to fix notations. Hereafter 
h=l. 

Given a symplectic vector space {V^ lo), where V has even dimension 2n and 
w is a nondegenerate skew-symmetric bilinear form on it, a Weyl system is a 
strongly continuous map from V to the group of unitary operators on some 
Hilbert space H: 

W : V ^U{H) (13) 

satisfying the condition 

W{vi)W{v2)W\vi + V2) = le^'^^"^-''^^ . (14) 

It is a projective unitary representation of the Abelian vector group associated 
with V. 

A theorem due to von Neumann [35] establishes that such a map exists for 
any finite dimensional symplectic vector space. Indeed, the Hilbert space Ti 
can be realized as the space of square integrable functions on any Lagrangian 
subspace of V. By using a Lagrangian subspace L, this is a subspace L of 
dimension half of the dimension of the space V and such that u){x,y) = for 
all x,y & L, and its dual L* it is possible to decompose V into V ^ L (B L* = 
T*L. Because L is Lagrangian, Eq. implies that the unitary operators 

corresponding to these vectors will commute. We consider the Lebesgue measure 
on L and we construct the following specific realization of W. The action of the 
restrictions W^Ilj W^Il* of the map W to the subspaces L and L* , on C^{L, cTx) 
is given by: 

{W\L{yWix)=H^ + y), (M^|L.(a)^)(x)=e^"(^V(x), (15) 

for x,y e L, a & L*, ip G C'^{L,d"x). 

The strong continuity requirement in the definition of W allows to use Stone 
theorem to get 

p^(to) = 6**^^") , yveV, (16) 

with R{v) the infinitesimal generator of the one parameter unitary group W{tv), 
t G M, depending linearly on v. 

In the following we shall be mainly concerned with the simple case V — M.'^, 
where, for simplicity, we introduce {q,p) coordinates. Then, by denoting 

P = R{1,0); Q = i?(0,l) (17) 
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we may also represent the Weyl map in the famihar form v = {q,p) — > W{v) = 
D{q,p), where the displacement operator D(q,p) is given by the formula: 

Diq,p)^e^p[i{pQ-qP)], (18) 

which turns out to be an irreducible representation of the Weyl-Heisenberg 
group (see section [5751 for more details on this topic). 

By introducing complex coordinates z ~ -^(q + ip) for V, we may write the 
displacement operator, eq. (|18p . in terms of creation and annihilation operators 
as: 

D{z) = exp (za^ — z*aj . (19) 

Generalizations of previous formulae to the case V — M^" are obvious, just 
by considering n— dimensional vector operators. 

Thus, a Weyl system provides, for any choice of a Lagrangian subspace, 
canonical pairs of operators and the displacements operators, i. e. a projective 
irreducible representation of the translation group. It is now possible to asso- 
ciate an operator with any function / on V admitting a (symplectic) Fourier 
transform /. Consider 

f{q,p)^ J (r'ad"xf{a,x)exp[i{aq- xp)]. (20) 

The operator W{f) associated to / is obtained by a clever use of the Fourier 
transform, substituting the (symplectic) Fourier kernel exTp[i{aq — xp)] with the 
Weyl displacement operator D(a, x) = exp [i [aQ — xP)] : 

W{f) = j d"ad"xf{a, x) exp [i {aQ - xP)] . (21) 

The above formula defines a unitary isomorphism between the space of the 
Hilbert-Schmidt operators on C^{L, d"x) and the square integrable functions of 
/;2(L©L*,d"a;d"a). 



3.2 Wigner functions 

We have considered a Weyl system to be a projective unitary representation of 
an Abelian vector group V of even dimension. Another useful interpretation of 
a Weyl system comes from the following considerations. 

Consider a fiducial vector IiJjq) in the Hilbert space Ti. carrying the projective 
unitary representation of V. We may consider an immersion of V into Ti. by 
means of the map 

V 3v< — >\v)^ W{v) IV^o) ■ (22) 

We denote this map by Wo : V H. The image of Wo is a submanifold of H, 
it is not a subspace. If we consider a fiducial operator Aq £ End{H) , we may 
in a similar way immerse V into the space of operators acting on Ti. by setting 

Wo : V" ^ End(H) , v i — > W\v)AoW{v) =: A{v) . 
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As it happens with any immersion of a manifold M into a manifold TV, we can 
consider the pull-back to V of the covariant tensor fields on H. As a matter of 
fact, on account of the probabilistic interpretation of quantum mechanics, the 
immersion of V should be considered to take place into the manifold of rays, 
i.e. TZ{Ti) the manifold of rays in Ti. or the complex projective space associated 
with 7i. 

With any immersion (j) : Ai — > A/" we have a map 0* : T [N) ^i-M.) 
which, however, need not be surjective. Moreover, dealing with infinite di- 
mensional manifolds {Ti or End(7i)) the properties of the pulled-back tensors 
depend on the specific immersion we deal with. In particular we may require 
the map v i — > \v) , or v i — A{v), to satisfy appropriate measurability, con- 
tinuity or differentiability properties. Depending on the use of the pulled-back 
tensors we may prefer one immersion over the other and, for instance, prefer 
the Schrodinger picture over the Heisenberg picture. 

By using the decomposition of V into different Lagrangian subspaces, say 
L X L' (we use the cartesian product notation to stress we are considering it as 
a manifold rather than as a vector space), we may consider 'eigenvectors' of the 
position operators, say \q) or \q') . In this way we may 'pull-back' any vector \tjj) 
to a wave function ip^q) = {q\4') or ipil') = ('Z'lV') ! similarly for an observable 
A we have the 'matrix coefficient' {q \A\ q') ~ fA{Q, l')- 

Had we chosen the decomposition V = L x L* , i.e. using 'eigenvectors' \q) 
and \p), we would have had (q or {p\A\q) , classically this corresponds to 

the use of boundary values, q and q\ or initial Cauchy data, q and p. These 
various matrix coefficient functions are connected by means of the completeness 
relations 

J \q)dq{q\=I = J \p) dp {p\ . (23) 

This general idea of considering the representation as an immersion of a 
manifold V into the Hilbert space H [37] , allows to consider the pull-back of the 
algebra structure on the operators. We may define a star-product by setting 

(t>*{A) * (j)*{B) (t)*{AB) . (24) 

In this way, on the subspace of functions we obtain by means of the pull-back 
of the operators a nonlocal and noncommutative product. 
In the Schrodinger picture we would have 

fAiv) = ^ , (25) 
{v\v) 

we could also consider 

fA{v) = {i^o\W\v)AoW{v)\iJo) , (26) 

and so on. In general, these functions are called symbols of the corresponding 
operators and carry a specific qualification to keep track of the specific immer- 
sion one is considering (we have Weyl symbols, Berezin symbols, Hormander 
symbols, and so on). 
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At this point it should be clear that if we deal also with unbounded operators 
on Ti. we end up with a variety of situations already at the level of topologies 
we are willing to consider on End(7i). Additional problems will arise from the 
specific choice of the fiducial vector we start with (if we need to take derivatives, 
we better deal with smooth or analytic vectors ^Sj). Usually, it is better to 
leave some of these choices unsettled and take them into full account only in 
each specific problem. 

Having clarified some aspects connected with the pull-back of observables, 
let us turn now to the pull-back of states. We recall that states are usually 
considered to be normalized positive functionals on the space of observables. 
States are not a vector space but we may consider convex combinations. Pure 
states are those that cannot be written as convex combination. To avoid some 
pathologies, very often states are also required to be normal. 

Thanks to Gleason's theorem, states are also called density operators, how- 
ever this may be misleading because it may give the impression that they should 
be considered with the same topologies as the operator algebra. However the 
star-product we have considered allows to distinguish the algebra, associ- 
ated with operators, which acts on the space of states. Thus, while the 
pull-back of states or observables always provides us with functions on V, the 
subsets to which they belong enjoy quite different properties and therefore it is 
advisable to avoid considering them as mathematical entities of the same kind. 

The distinction will reappear to be crucial when we would like to compare 
with the analogous situation in classical mechanics described on the same phase 
space V. Again here states and observables are quite different, pure classical 
states end up being distributions, while observables, due to the possibility of 
taking Poisson brackets, are usually required to be smooth functions. This 
distinction plays a very relevant role when we consider their associated Radon 
transforms to represent tomographic states or tomographic observables. 

We may now define a Wigner function the way it was considered by Wigner. 
Given a state {ip) we form the rank-one projector 

which defines a normalized positive functional on the space of observables. We 
may consider, as we have stressed, either the density state 

or the function 

i.e. the matrix coefficients of the rank-one projector in the (position, position) 
representation or in the (position, momentum) representation. By using the 
completeness relation we have 

i^W)^ I {i^\p)dp{p\q') (30) 
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which provides us with the transformation from one representation to another. 

It is now feasible to use convex combinations of pure states to define generic 
states and their associated Wigner representations. The Wigner function yV(p, q) 
of a density state p{q, q') is defined, restoring the Planck constant /i, as: 



and it results 



W{q,p) = 2exp(-2p(7) / W(g',p')exp 



{q'p' - 2pq' - 2qp') 



dq'dp' 



(31) 



(32) 



Now a simple manipulation, changing the variable x/2 — s and extracting from 
the integral the bra {'ip\ and ket ji/'), we get: 



W((7,p) 



j e-2»P^/'^|g-s)(g + s|ds 



1^) 



(33) 



3.3 Transformation properties of Wigner functions 

We would like to understand the transformation properties of Wigner functions 
under linear symplectic maps and dilations on phase space. The exploration of 
these issues will suggest the possibility of extending the definition of Wigner's 
function to the space of irreducible representations of the Weyl-Heisenberg 
group and to discern their homogeneity dependence on Planck's constant h. 
We will make these comments precise in what follows. 

The Weyl map allows to associate automorphisms i/^ on the space of uni- 
tary operators with elements of the symplectic linear group Sp{V,Ljj) of V, 
according to the following diagram 



V 



V 



w 



w 



uin) 



(34) 



for any e Sp{V,uj), by setting 

= W{^{v)) = ulw{v)U^, Vw e y . (35) 

Recall that an isomorphism <j):V is symplectic if w((^(u), 0(w)) = uj{u,v) 
for all u, w e Moreover, the group of all symplectic isomorphisms of V can be 
identified with the matrix symplectic group Sp{n) by choosing a symplectic basis 
on V . In other words, the automorphism v^^ corresponding to the symplectic 
linear transformation (j) of V, is a inner automorphism of the group of unitary 
operators, that is there exists a unitary operator Uip such that I'tpiV) — U^VU^ 
for all V e U{Ti)^ because it belongs to the connected component of the identity 
of the automorphism group [39] . At the level of the infinitesimal generators of 
the unitary group, we have 



UlR{v)U^ = Ricbiv)) 



(36) 
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Further insight on the physical meaning of the Wigner function of a density 
state p, was obtained from its representation as the expectation value of the 
shifted parity operator P{q,p) (see, e.g., Royer [IS])- In fact, we can write the 
expression given by Eq. (|33p of the Wigner function yV{q,p) corresponding to 
the state \^), as: 

W(<z,p) = 2(V|P(g,p)|V), (37) 
with V^qjp) being the shifted parity operator: 

r{q,p)= J e-''^P^/^\q-s){q + s\ds. (38) 

Notice that 7^(0, 0) = /|— s)(s|(is is just the parity operator defined as (7'?/') (g) = 
ip{—q) or, equivalently the unitary operator V satisfying: 

VQV = ~Q, VPV = -P. (39) 

Now, we get immediately that: 

V{q,p)^D{q,p)VD{q,p)^, (40) 

where the displacement operators D{q,p) have the usual form given by Eq. ([T5)) . 
Then, the Wigner function corresponding to the pure state \ip) can be readily 
written in the form: 

Wiq,p) = 2 (V^l D{q,p)VDiq,p)^ |^) , (41) 

and for a given density state p we obtain: 

Wpip, q) - 2Tr [pDip, q)VD\p, q)] = 2Tr [pD{2p, 2q)V] . (42) 

Because the displacement operators provide a specific irreducible represen- 
tation of the Weyl-Heisenberg group, the previous formula makes apparent the 
possibility of generalizing Wigner 's function as a function on the space of ir- 
reducible representations of the Weyl-Heisenberg group. In fact the Weyl- 
Heisenberg group WH(n), for n = 1, may be presented as the group of triples 
of real numbers (p, q, t) with the composition law: 

{q,p,t)o{q',p',t')^{q + q',p + p',t + t' + ^ {pq' - qp')). (43) 

The associated canonical operators with their commutation relations [Q, P] = 
il are a realization of the Lie algebra of the Weyl-Heisenberg group, and a 
irreducible unitary representation is provided by: 

U{q,p,t)^D{q,p)^". (44) 

In the general case n > 1, the irreducible representations of WH(n) are 
parametrized up to a unitary equivalence by a real parameter 7 [41] . Kirillov's 
theory of coadjoint orbits [42j provides a natural way to construct them. In fact. 
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Kirillov's theorem establishes that for nilpotent groups there is a one-to-one 
correspondence between coadjoint orbits of the group and equivalence classes of 
unitary irreducible representations of it. It is easy to check that for the Weyl- 
Heisenberg group the space of coadjoint orbits has two strata, the regular one 
whose coadjoint orbits are copies of the symplectic linear space {V, ui) and are 
labelled by 7 ^ 0, and the singular stratum, corresponding to the label 7 = 
whose coadjoint orbits are points, hence giving rise to trivial representations. 
The parameter 7 weights the central element of the group and it can be easily 
read out from a given irreducible representation looking at U^{0,t) — e''''* and 
therefore the action of Aut(WH(n)) on the set of irreducible representations 
can be analyzed. 

In order to introduce a generalized notion of Wigner functions [43] for repre- 
sentations with 7 7^ 1, we have to choose previously a representative Uj out of 
any equivalence class [U]^ . We choose, for 71 = 1, the representatives for 7 > 
as: 

U^{q,P,t) C/^=i(v^g, V7P,70 = D{^q, ^q)e'-'K (45) 

Once a representation Uj has been chosen, the parity operator V given by 
Eq. p9|) may be expressed as: 

P.2/^„(^„^rt.i/M:z„„rt, ,46, 

^From this expression the properties: 

V = V^; VU^iq,P,t)V = U^{-q,-p,t), (47) 

readily follow. 

Now, given 7, we define the associated (generalized) Wigner function of a 
density state p as 

Wpiq,pn) : =2TT[pU^iq,p,t)rUl{q,p,t)] 
= 2TT[pD{2^q,2^p)V] 

= Wpi^q,^p;l). (48) 

We remark that, while the dependence on the parameter t disappears and the 
function is invariant on the subgroup (0,t), a new dependence on the represen- 
tation label 7 appears. 

We now consider the action of a dilation (l)x : 0a (9,^,0 = {M^ ^Pt ^^t)- 
Then, as a result of our choice of the representatives U^, we get 

[/^(0A(g,P,t)) = U^{Xq,XpA''t) = Ux2^{q,p,t). (49) 

So, the Wigner function transforms as: 

Wp(0A(g,p);7) = Wp{Xq,Xp;j) 

- 2TT[pUx2^{q,p,t)VUx2^{q,p,t)] 

= Wp{q,p;X^j), (50) 
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while 



X^^dqdp 
2^ 



Wp(Ag,Ap;7) 



X^jdqdp 
2^ 



Wp{q,p;X^-f)^TTp. (51) 



The dilation transformation may be more interestingly written as: 

Wp[Xq,Xp;^) =Wpiq,p;j). (52) 

We observe that the dilation {Xq, Xp, X^t) yields the expected dilation on 
the label 7, which is 'dual' of the parameter t. For an infinitesimal dilation 
A = 1 + e we may expand: 



Wp{q,P-n) = Wp^l + e){q,p)- 
= Wp(g,p;7) + e 



7 



(1 + ^) 
9W, 



(53) 



P, , dWp ■ 
^^-(9,P;7)-27^(g,p;7) 



0(e^ 



where we have used the notation v = {q,p) and 



dWp dWp dWn 
— q—^ +p- 



9v dq dp 
Then we obtain the following differential equation for the Wigner function: 



(54) 



9W, 



^(v;7)-27^(v;7) = 0. 



(55) 



So far, we have put h = 1. It is possible however to study the dependence 
on h by using the displacement operators given, instead of Eq. (fT8| . by the 
expressions: 

and the canonical commutation relations: 



Dh{q,p) = exp 



(56) 



1 



[Q,P] = z/, 



(57) 



while t gives place to t/h and the unitary representation given by Eq. (|44p 
becomes: 



(58) 



so that eventually we get the above formulae with 7 replaced by "f /H everywhere. 
In particular, for the Wigner function we have: 



Wp(V^v;l)=Wp(v;7) 



Wo 



v;l 



(59) 
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Under the action of a dilation, 

(Av, Xh) ^ (Av; 7/A2) ^ (Av; j/XHi) (60) 

and we may choose 7 = 1, to get a differential equation for the Wigner function 
Wp (v; i) corresponding to the infinitesimal 'dilation' (Av; 1/X^fi). 

Notice that the scaling properties Eq.s ([551 ISO)) is consistent with the de- 
pendence on h of the Wigner function. We recall that the density state p has 
the dimension of an inverse length , where £ ~ s/hj^/rmj : so p{x,x') = 
i~^p'{x/£,x' /£). Then it is easy to check that the following property holds: 

We refer to [43] for more details and further results on this direction. 

3.4 Tomograms, Wigner functions and Radon Transform 

By means of the reconstruction formula ([T2]) the Wigner function W{p, q) of a 
density state p{q, q') may be recast in the form: 

W{p, 9) = — / T{X, p, v) exp [i{X - pq- vp)] dXdpdu. (62) 
27r J 

The above equation explicitely contains the Planck constant to be coherent 
with Eq. ([T2| . Hereafter h= 1, however. We recall that in general the Wigner 
function is not a fair probability distribution as it is not non-negative; never- 
theless, it is a function on the phase space of the system, and its reconstruction 
formula (|62p is just the Radon anti-transform of the tomogram. 

The Radon transform [10] originally was formulated to solve the problem 
of reconstructing a function /(p, q) from its integrals on arbitrary straight lines 
pq -\- vp = X 'm the {q,p)-^\a.ne 

J f{p,q)5{X - pq- vp)dpdq=: {nf){X,p,y). (63) 

Here 5 is the Dirac delta function and the parameters X, p, v are real. The 
homogeneity property follows from the properties of the delta function. The 
inverse transform reads: 

/(p, q) = / (7^/)(X, p, y) exp [i(X - pq - vp)] dXdpdu. (64) 

(27r) J 

Remark: Additional hypotheses, such as global integrability conditions, are 
required to guarantee the uniqueness of the inverse transform [33] . We point out 
that here a subclass of functions is selected by requiring that our 'manipulations' 
provide us with an injective map. 

In a general sense, we may call the Radon transform (|63p . the tomogram 
of the function f{p,q), therefore the ambient space for tomographic states is 
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provided by the range of the Radon transform when it is apphed to properly 
chosen functions on phase space. 

The problem we address now is the following. Let us select two classes of 
functions on the phase space satisfying special conditions. The first class of 
functions consists of all probability distribution densities on phase space (the 
q — p plane) describing states of classical particles. The second class consists 
of all the Wigner functions describing quantum states thought of as rank-one 
projectors. We study tomogram properties of these two classes. There exists 
also a class of tomograms which are Radon transforms of the Weyl symbol of 
observables. These tomograms are not probability densities. 

So, the symplectic tomogram we dealt with in section [5] may be eventually 
interpreted as the Radon transform of the Wigner function 

T{X,tx,y) = j l-W{p^q)5{X-fiq-vp)dpdq (65) 

= j p{yiV')'^*x^viy)vxt,v{y')dydy' . 

The standard description of classical states with fluctuations is given by a 
non-negative joint probability distribution function /(p, q) on the phase space 
(a plane, for a particle with one degree of freedom). The function is normalized, 
i.e. 

j f{p,q)dpdq^l . (66) 
The classical state tomogram {TZf){X, /i, can be written in the form 

inf){x,^i,ly) = {^ix-^Iq-l^p))^ (67) 

where the average is done using the probability distribution /(p, q) in the phase 
space |311 132) . The tomogram is the probability distribution function in a ro- 
tated and scaled reference frame on the phase space. It can be expressed in 
terms of a scaling parameter s and a rotation parameter 9 : 

fi — scos9, v~s^^sm9. (68) 

For fixed fi and v one then gets a line X = fiq + vp in the plane {q,p) with an 
orientation 9 from the position axis. Thus the physical meaning of the variable 
X is that it is the 'position' of the particle measured in the reference frame of 
the phase-space whose axes are rotated by an angle 9 with respect to the old 
reference frame, after preliminary canonical scaling of the initial position q ^ sq 
and momentum p —> s~^p. The coordinates X and Y = —s^h'q + s~^^p provide 
a canonical transformation preserving the symplectic form in the phase space. 
For that reason the classical tomogram is called 'symplectic'. 

In the quantum case, Eq. I|65p can be written in a form similar to Eq. (|67p : 

TpiX,^I,l^) = {S{X-^lQ-,yP))^. (69) 
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The difference with Eq. ([S7| is that here the position and momentum are quan- 
tum operators Q and P, and therefore we have to take into account uncertainty 
relations. The averaging in Eq. is done using a density state p, i.e. 



(70) 



For fixed jj, and the operator 



(71) 



together with its conjugate 



(72) 



satisfies the canonical commutation relations of the Weyl-Heisenberg algebra: 
[X{^,,v),Y{^,,v)] = [Q,P]. The observable X{fj,,i') is a new position operator, 
i.e. the position after a symplectic (linear canonical) transformation in the quan- 
tum non-commutative phase space {Q, P) of the particle. The real variable X 
gives the possible results of a measure of X (/i, i^) and runs over the spectrum of 
X (n, u) . In this way a description of quantum tomograms is recovered in com- 
plete analogy with the classical case. So the tomogram is also 'symplectic' in the 
quantum case, it is associated with an automorphism of the Weyl-Heisenberg 
algebra. 

In classical mechanics, the transition from the distribution function of two 
canonically conjugate variables (position q and momentum p) to the distribu- 
tion function of a (/i, ;/)— family of position variables (position X) does not play 
a crucial role, due to absence of quantum mechanical constraints like the uncer- 
tainty relations of Heisenberg [JS] and Schrodinger- Robertson Wf\ SHI SH] • 
On the contrary, the use of tomograms in quantum mechanics should encode 
the properties required to allow for the uncertainty relations. 

So, the tomograms of all admissible functions /(p, q) form an ambient space 
(here admissibility means only that the Radon transform exists and is one- 
to-one). This space contains the subset of probability densities. In turn, the 
subset of the probability densities contains two subsets. One subset contains the 
Radon transforms of Wigner functions which are probability densities (quantum 
domain). The other one contains the Radon transforms of classical probability 
distributions on phase space. These two subsets have a not empty intersection. 
Both these subsets are embedded into the total set of tomograms, which there- 
fore contains tomographic functions corresponding neither to classical nor to 
quantum states. 

3.5 Distributions and quasi-distributions: classical and 
quantum 

As we have argued, on the same space an object like a symplectic tomogram 
T(X, /i, v) may determine a state both in classical and in quantum domain. Let 
us discuss some difference which exist for these two domains in the context of 
the tomographic description. State tomograms in both domains must satisfy 
the following common requirements: 
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1. Nonnegativity: T(X, /i, > 0. 

2. Integrability: J T{X, ^,iy)dX < oo. 

3. Homogeneity: T(AX, A/i, Ai^) = jj|-T(X, /i, i^). 

Other properties of the symplcctic state tomograms are required to distin- 
guish the states in quantum and classical domains. For example, the necessary 
condition for the tomogram of a classical state is the nonnegativity of its Radon 
anti-transform, i. e. 

J T (X, ^, v) exp \i{X - M - vv)\ dXd^idv > . (73) 

The violation of this inequality means that the tomogram does not describe a 
classical state. On the other side the condition for a symplectic tomogram to 
describe a quantum state can be formally written in an analogous way, as 



J T{X, fi, u) exp [i{X ~ nQ- i^P)] dXdfidiy > 



(74) 



which means that the operator obtained from the above Radon anti-transform, 
being a density state, must be a non-negative normalized functional on observ- 
ables. If the inequality is violated, the tomogram does not describe a quantum 
state. Tomograms satisfying both conditions belong to the intersection of quan- 
tum and classical domain, so they may be chosen as starting Cauchy datum of 
either quantum or classical time evolutions. For example, the gaussian tomo- 
gram of coherent or squeezed and correlated states [121 [ST] satisfies both 
inequalities. 

Now, the question arises if the previous requirements of nonnegativity, in- 
tegrability and homogeneity are sufficient to select only classical or quantum 
tomograms. In other words, there exist tomograms satisfying the three require- 
ments, but violating both inequalities? It would mean that these tomograms 
need not describe a state, neither a classical one nor a quantum one if additional 
requirements are not met. An example of such a kind of tomograms may be 
manufactured to answer in the affirmative this question. The example is pro- 
vided by scaling the parameters /x, ly of the tomogram of the first excited state 
of an harmonic oscillator 



2e 'i^^T^ 

Ti (X, /i, y) = -^—^ (75) 



by means of a real parameter A obtaining 



2e (ITo^Tp^^'' X"^ 
Tx {X, , .2 .2 ■ (76) 



This new tomogram is still positive, integrable and homogeneous. But the quan- 
tum inequality is not fulfilled by Tx for A 7^ 1, as discussed in [52] where, in a 
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different contest, it is shown that the fidehty of the scaled tomogram and the gen- 
uine tomogram of the harmonic oscillator ground state is negative: / = —2 |A|^ 
for small A. Also, the classical inequality is violated by the scaled tomogram, 
because its Radon anti-transform yields a distribution (generalized function) 
which is negative for small g, p. 

In conclusion, the set of non-negative, integrable and homogeneous tomo- 
graphic functions is divided into three parts: one containing tomograms of quan- 
tum states, another one containing those of classical states and a third part 
containing neither. The first two parts intersect each other in the domain of to- 
mograms satisfying both quantum and classical inequalities. The third part does 
not intersect the others, and contains tomographic functions describing neither 
quantum nor classical states therefore they violate both inequalities. As a con- 
sequence, more constraints than the previous ones are needed to unambiguously 
select quantum, or classical, tomograms in order to give a tomographic version 
of quantum mechanics fully equivalent to the usual formulations and yielding 
the classical mechanics in an appropriate limit. At least in principle, for the 
symplectic case they are sufficiently described by the quantum, or classical, 
inequality. 

As intrinsic characterization of classical tomogram is the property that the 
Fourier components of the classical tomograms is non- negative (see Eq. (1751) ). 
The property of quantum tomograms given by inequality ([71]) is also intrinsic 
but needs the additional construction of operators Q and P on a Hilbert space. 
One can formulate the inequality in a form equivalent to Eq. (|74[) without 
introducing these operators. This is accomplished replacing the nonnegativity 
condition of the operator with Sylvester criterion for the matrix of the operator 
corresponding to the integral kernel of Eq. ([74|) . It means that in any basis 
the principal minors of such a matrix are non-negative. This condition is given 
in form of algebraic integral inequalities, not containing formally mention of a 
Hilbert space. 

Another intrinsic description of the set of quantum tomograms can be for- 
mulated in the following way. The tomograms of coherent states are gaussian 
and belong to the intersection of quantum and classical sets. The intersection 
is invariant under the operation of taking convex sums of the tomograms. Due 
to the completeness property of coherent states, the composition of all coherent 
state tomograms given by Eq.s ()137p . ()138p provides the set of all quantum to- 
mograms corresponding to pure quantum states. Then, combining with convex 
sums these pure quantum states, we get any mixed state tomogram. In this way 
the whole set of quantum tomograms is recovered. 

3.6 Tomographic sets and state reconstruction 

^From the conceptual point of view, Pauli's problem raises a central question in 
the formulation of quantum mechanics. In general, any system may be described 
by a set O of observables and a 'dual' set of density states S, which together 
give rise to probability measures on the real axis, they are just the probability 
distributions of the values of the observables in the states. 
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In the usual formulation of quantum mechanics, any observable A, i.e., an 
hermitian operator, uniquely determines a projector valued measure (PVM) 
Pa{E) on the sets E of the Borel cr-algebra of the real line [521 ISl], so that from 
a density state p a probability measure mp,A can be defined as 

mp^A{E):^TY{pPA{E)) . (77) 

It is just the probability that the mean value of the observable A in the state p 
belongs to E. As a consequence: 

mp,^(R) - 1 . (78) 

The mean value of A on the state p, when it exists, can be written as an integral 
over a real variable A with respect to that probability measure: 

Tr(pA) - j Xnip^A ■ (79) 

This for pure states p = \ip) reads 

(^|A|V> = f Xmp,A . (80) 



and a functional calculus for hermitian operators can be constructed by defining 
the operator f{A) as: 

{^P\fiA)\yj) ^ I f{X)mp,A . (81) 

for any integrable function /. So, the knowledge of mp^A for any state p and 
fixed A allows for the reconstruction of {ip \A\ ip) for all ip in the domain of A 
and therefore of A. In fact, the matrix elements of A in any chosen basis of its 
domain are given by the polarization identity: 



{^P\ip)^i \\^j + ^\f + U-^\f-tU + iip\f+iU-i^\f . (82) 



Viceversa, when p is fixed, the knowledge of nip, a for any observable A, in 
particular for all projectors , allows for the reconstruction of Tr(p - 

{ijj \p\ ijj) and therefore, by polarization, of p. 

For instance, when A is the position Q, the associated projectors Pq{E) act 
on wave functions as a multiplication by the characteristic function of the Borel 
set E: 

(Pq{E)^){x) ^ Xe{x)^{x) (83) 

Now, a density state p can be spectrally decomposed in terms of rank-one pro- 
jectors as (because selfadjoint compact operator): 

p = X!fc ^^^^ ' - ^ ^' (^'^) 
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therefore the previous formulae became: 

mp,Q{E) := Tr{pPQ{E)) = J2,^kTr{PkPQ{E)) 

= ^^oik J XE{x)\ipkix)\'^ dx ^^^ak j \ipk{x)fdx 

= j ^^Qfc |V'fc(a;)|^rfx = j p{x,x)dx. (85) 

This shows that the probabihty measure mp^Q{E) is absolutely continuous with 
respect to the Lebesgue measure on M, with Radon-Nykodim derivative 

p{x, x) = ak \il>k {x) I ^ (86) 

which is just the diagonal part of the density matrix in the position representa- 
tion. Then the mean value of the position operator is: 

Tr(pQ) = Jxp{x,x)dx (87) 

or, for pure states \tp): 

(tplQl^p) = I x\^{x)\''dx . (88) 



Thus, Pauli's problem may be reformulated as: To determine the state from 
the knowledge of a pair of probability measures rUp^A, i-e. when A is the po- 
sition Q or the momentum P. This set is not sufficient, while the symplectic 
tomography provides a set of obscrvablcs X(/i, v) — fiQ + vP which is sufficient 
for the reconstruction. Notably, the symplectic set is generated by the posi- 
tion operator Q under the action of a family of unitary transformations S{iJi, v). 
Introducing the auxiliary parameters A, 9 as 

= e-^cos6' = e~-^sin^ , (89) 

we can write 



S{p,, y) = exp 

so that 



^(QP + PQ) 



exp 



(90) 



S{ix, v)QS^ in, v)=p.Q + vP . (91) 

Thus, the transformation ^(/x, v) yields the appropriate probability measure 
associated to the observable X{pb,v) 

rUp^xMiE) := Tv{pS{i,,u)Pq{E)SHi,,u)) = f {X \S\i^,i^)pS{i^,i^)\ X) dX 

JE 

(92) 

whose density is just the tomographic probability distribution 

Tp{X,^l,l^) -.^ {X\S^{ix,v)pS{ii,v)\X) = Tr(pS{^Ji,v)\X) {X\S\ii,v)) (93) 
where the kets \X) are the eigenkets of Q : Q \X) = X \X) . 
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3.7 Rank-one projectors as tomographic sets 

In view of formula (|93p we may consider in general tomographic, i.e. (possibly 
over-) complete, sets of rank-one projectors. In a sense, they are the elementary 
'building blocks' of any tomography. Moreover, as we will show in the following, 
tomographic family of rank-one projectors allow to clarify readily the ingredients 
of a tomographic reconstruction formula. 

We start with an abstract finite dimensional case. Assume the Hilbert space 
of the vector states Ti to be n— dimensional, so that rank-one projectors span 
a n^— dimensional Hilbert space H =7i ® Ti,* , containing all the density states 
as well as the (bounded) operators on Ti., i.e. H ~ B{T-L). The scalar product 
is given by the trace: {A\B) — Tr{A^ B). A minimal tomographic set is a basis 
{Pfc} , fc € {l, n^}, of rank-one projectors, which may be orthonormalized by 
a Gram-Schmidt procedure: 



k=l 

In general, every \Vj) is a linear combination of projectors, rather than a single 
projector like \Pk) ■ Then a resolution of the super-unity on H in terms of the 
P's reads as 



El^')(^«l= E 7*,7.,-P,Tr(Pz.) 



1=1 i,j,i=i 

2 2 

n n 

1=1 J=l 

where the dual set of Gram-Schmidt operators {Gk} ,k£ has been 

introduced: 

2 2 

n n 

\Gi)^J2^:im= (95) 

i=l i,j = l 

We observe that G; is a nonlinear function of the projectors {Pk}, because also 
the coefficients 7's depend on the projectors. Moreover: 

2 2 

71 71 

j=i ],k=i 

2 

n 

= (96) 

Thus, Eq. ([94]) shows that a resolution of the (super-) identity, associated 
with a tomographic reconstruction formula, is determined by a pair of dual sets, 
the P's and the G's. Besides the role of the dual operators may be interchanged 
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in such a formula, so that a tomography would be better defined in terms of a 
pair of dual sets. In the context of the harmonic analysis and of the wavelet 
signal analysis, these dual sets are known as dual frames [55] , 

Formulae similar to the previous ones hold for any tomographic, i.e. (over-) 
complete, set of rank-one projectors. In fact, along with the minimal tomo- 
graphic set discussed above, that is a quorum of rank-one projectors, it is very 
useful to deal with over-complete or even maximal sets of rank-one projectors, 
obtained for instance by acting on a fiducial projector Pq with a unitary repre- 
sentation of a Lie group. The previous summations became then integrations 
over the orbit fl through Pq. 

A suitable illustration of such tomographic sets is given by the qubit (spin 
1/2) tomography over a 2— dimensional Hilbert space of the vector states H. Out 
of the standard basis vectors where m = ±1/2 is the spin projection on 

the z-xis, a fiducial projector Pq = \m){m\ is rotated by means of the operators 
U of an irreducible representation of the group SU (2) : the tomogram of any 
operator p is defined as 

Tp{m,U) = Tr([/|m)(7n|[/V) = {m\U'' pU\m) . 

Then fl is the orbit of the co-adjoint group action in the dual of the algebra, that 
is the Bloch sphere of all rank-one projectors. The projector U\m){m\U^ , 
with U parametrized by the usual Euler angles {6, ip, -0), corresponds to a point 
on S^determined by a unit vector n = (sin cos 0, cos cos </>, cos 0). In other 
words, we can use the (9, (p) parametrization to write a generic projector in 
matrix form as (m = 1/2) 



l-f cos 61 
e'"^ sin 6* 



e"*"^ sin( 
1 — cos 6 



(97) 



Then, the corresponding resolution of the identity reads [16] : 

|G(0,</>))(P(0,</>)|sin 



(98) 



where in matrix form: 



G(^,0) = - 



1-1- 3 cos 6* Se^'-Psin^ 
3e'"^ sin 9 1 - 3 cos ( 



(99) 



SO that, for any density state p we get the dual reconstruction formulae: 

C2-K /'TT 



P = 



G{9, (f>)TT{P{9, (j))p) sin9d9d(l) 



Jo 

2-K I'-r; 



"'0 



P{9, <?!))Tr(G(6l, 0)p) sm9d9d(j) 



(100) 
(101) 



We remark that, however, in the infinite dimensional case the relation H = 
B{TL) is no more valid and there are several relevant spaces. In particular 
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the Hilbert space at our disposal is the space of the Hilbert-Schmidt operators 
H C B{H), which is the typical setting of the frame theory [56j, while the finest 
tomographic sets of rank-one projectors have to be complete both in such a 
Hilbert space and in the Banach space of the trace class operators. We will not 
insist here on the topological subtleties of the infinite dimensional case, they are 
discussed, e.g., in pTl [57] . 

3.8 General aspects of tomography 

The Pauli's problem and a possible solution have been previously analyzed 
within the machinery of spectral analysis of selfadjoint operators. In this sec- 
tion, inspired by that analysis, we provide a general setting for tomography 
together with some general considerations. 

When S is the set of states of a physical system and A a suitable subset of 
the observables O, a tomography T is a map from Sx [A C O) into the set of 
probability measures on the real line M. It is required that T is such that if the 
probability measures A) are known for all A g it is possible, at least in 
principle, to reconstruct p. 

If ^ = O a tomography is available, via spectral analysis, as shown previ- 
ously. In this case A is a huge linear space of selfadjoint operators, bounded or 
not. 

However, as we have seen, a tomography is available even if A is restricted 
to a subset C C O of all the rank-one projectors. Now O' is a unit spherical 
surface in the Hilbert space, in general infinite dimensional. It may be more 
useful to restrict ^ to a subset of O specified by some (multi-) parameter fj, 
which varies in some index set M . A must still be such that the reconstruction 
of any state p is possible. The tomogram T{p, A) appears now as T{p, p){E) = 
T{p,Af,){E) = Tt{pPa^{E)) and is a probability measure: < T{p,p.){E) < 1, 
for any Borel set E of reals. 

In principle the elements in A can be selfadjoint operators with different 
spectra but it is more convenient to deal with iso-spectral operators, so that 
a spectrum cr C M is associated with A. In this case, for any p, p, we have 
that T{p, p){E) = for all sets E which do not intersect the spectrum cr and 
T{p, p){E) — 1 for sets E containing the spectrum a. 

If (T is a purely continuous spectrum, as in the symplectic tomography, and 
T{p, p) is absolutely continuous with respect to the Lebesgue measure, T{p, p) 
will appear as p, p)dX, where X e cr is the spectral variable, with the 

property T{X, p, p)dX = 1. In other terms, T{X, p, p) is a probability density 
function on the real line concentrated on the spectrum cr. 

If CT is a pure point spectrum, T{p, p) will be concentrated on the points Xk of 
cr, with T(p, p){Xk) = Tv{pPx^.fj) — T{k, p, p) a positive function of k such that 

T(fc, p, p) — 1. Now {T(fc, p, p)}j, can be regarded as a probability vector T 
with a finite or countable number of components, i.e. a vector with non-negative 
components whose sum is one. When the spectrum is non-degenerate finite, 
the number n of components of T is just the dimension of the Hilbert space. 
Geometrically, these n components may be considered as the coordinates {rfc} 
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of a point belonging to the simplex ti + T2 + • • • + t„ = 1 of R". If the spectrum 
is degenerate, the number of components of T is less than the dimension of the 
Hilbert space. For instance, if A is composed by rank-one projectors, a = {0, 1} 
and T has two components. However, we notice that the qubit tomographic 
set presented in the previous subsection is a non-degenerate case: given the 
component ti = Tp{m, u), the other one is T2 = \ — Tp{m, u) and the probability 
vector belongs to the simplex ri + T2 = 1. 



3.9 Weyl systems and tomography 

As we have seen, tomographic families of iso-spectral observables can be ob- 
tained by conjugation of a fiducial observable Aq by a parameterized set of 
unitary operators J7^, which may be the elements of a unitary (possibly square 
integrable [22]) representation of some Lie group Q and ^ will then be a point 
on a finite dimensional manifold. In this section the most common tomographic 
schemes are obtained in an unified approach, by using a particular representa- 
tion of a Lie group Q. 

We discussed in section 13.31 how the Weyl map allows to associate inner 
automorphisms on the space of unitary operators with symplectic isomorphisms 
of the symplectic linear space {V,lo). Moreover if ^ is a Lie group and T: Q ^ 
Sp[V,uj) is a linear representation of the elements of the group by symplectic 
maps, then we can associate an inner automorphism v^g of the group of unitary 
operators to any element g of the Lie group by using the analog of Eq. (j35p : 



^tAW{v)) = W{Tgv) = UlW{v)Ug 



(102) 



where now Ug denotes the unitary operator associated to Tg realizing the inner 
autormorphism vt^ ■ In other words the automorphism i>Tg , corresponding to the 
symplectic linear transformation Tg of V, is a inner automorphism of the unitary 
operators, as it belongs to the component of the automorphism group which is 
connected with the identity [39] ■ At the level of the infinitesimal generators of 
the unitary group, we have 



UlR{v)Ug = R{TgV) 



(103) 



The equation above allows to make contact with most used tomographic 
schemes. For instance, by acting with a rotation of an angle 6 followed by a 
squeezing 



we have 







cos f 
— sin ^ 



sm t 
cos ( 



cos 6 



sin 9 

e^^ cos ( 



(104) 



Then, by 
and 



UgR{v)Ul - RMv) 



e'*'cos0 — e ^ sin 



(105) 
(106) 



25 



we obtain 



UgQUl 

UgPUl 



cos ( 



e ^ sin t 



Q 
p 



cost 



sin6'P 
■ cos OP 



In this way, bearing in mind Eq.s 
graphic family of operators Xlpi,!/) 



(107) 

1}, we recover the symplectic tomo- 



UgQUl = e^cos( 



■ sin 6IP = nQ + vP = S{n, iy)QSHp, v) = X(p, v) 



(108) 

We get a general, unified approach by using the inhomogeneous symplectic 
linear group, which is a semi-direct product of the translation group and the 
linear symplectic group. In that case we have only projective representations, 
whose generators are /, Q, P, \{P'^ + Q'^), \{PQ + QP),\{P'^ - Q^). Then most 
common tomographic schemes may be obtained by acting with the group on its 
generators as fiducial operators ^o- 

For instance, by acting on Q, a = M, we get 

D{q,p)Sia,fi,iy)QS^ia,fi,iy)D^q,p) (109) 
= D{q,p)i^Q + vP)DHq,p) = i^Q + vP) + i^q - VP)I 



with 



£, ~ fJ- cosh a + v sinh a: rj = v cosh a + p, sinh a , 



(110) 



and where, restoring the explicit dependence of /i, v on the parameters A, 9, the 
transformation S{a, /i, v) <-!■ S{a, A, 9) reads: 



S{a, A, 9) — exp 



exp 



^^iPQ + QP) 



exp 



9^' 



Remark: Equivalently we could have acted on P, a 
same tomographic family: 



, obtaining again the 
(111) 



i?((Z,p)5(a, -A, 9 - ^)PS\a, -A, 9 - ^)D\q,p) 
= D[q,p){^Q + ilP)D\q,p) = {^Q + i^P) + {^q - i^p)I . 

We may recover from the above tomographic family of observables the usual 
symplectic tomographic scheme [HI |6^ . We do that by quotienting the whole 
inhomogeneous group with respect to its subgroup generated by the displace- 
ment operators times the hyperbolic rotations, and choosing a section 



(7 = G,p = 0, a = 0. 
Besides, by acting on the number operator 



1 t 



N, No, 



(112) 



(113) 
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we get the photon number tomography [5T1 [5^ [55] on a section 



a = 0, A = 0, 61 = 0, 



(114) 



while on the other section 



0,p = 0,a = 0, 



(115) 



the so called squeeze tomography [53] is recovered. 

In the coherent state tomography the fiducial operator is the projector on 
the vacuum state: a |0) = 0, 



Ao- |0)(0|, a -{0,1} 



(116) 



and the tomographic family of operators is generated by the displacement op- 
erators depending on a complex parameter z, T>{z) — exp [za^ — z*a] , eq (I19p . 
or equivalently, two real parameters. Thus one obtains a tomogram which is 
formally equivalent to the Husimi function, while one of the dual reconstruc- 
tion formulae is the Sudarshan's diagonal coherent state representation of an 
operator [57] . 

We can also obtain spin tomography. For this, we observe that unitary 
irreducible representations of SU{2) of dimension m > 2 are a subgroup of 
SU (m), which in turn is isomorphic to S0{2m) n Sp{2m, R), the intersection of 
2m— dimensional rotations and linear symplectic transformations. Now, for def- 
initeness, assume = with coordinates (a;i, a;2,Pi,P2) in the Weyl diagram, 
with T a 2— dimensional unitary irreducible representations of SU(2), then for 
any g G SU(2): 

V ^ u{n) 

Tg I I ut, (117) 

V u{n) 

The matrices Tg are real, orthogonal symplectic 4x4 matrices while the unitary 
operators Ug corresponding to the inner automorphism i^Tg are generated by the 
hermitian operators 



Jl 



\{PlQ2-P2Ql) 



(118) 



2) 



{p? + QI)-Up^ + Q 



= -[H{Pi,Qi)-H{P2,Q2)] 

which satisfy the commutation relations [Jh, Jk] — ichkiJi - The operators Qa, Pa,ct 
1,2, generate the displacement operator of the Weyl system. 

In fact, choosing the lagrangian manifold to be {(xi, a;2)} , or {(p, </?)} in polar 
coordinates, then Ji acts as ^(id/dip) on 7i = £^(R^) and its eigenfunctions are 
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?/'(p)e™'^ so that its spectrum is the set of the integer and semi-integer numbers. 
J3 is a semi- difference of two isotropic harmonic oscillators, so it has the same 
spectrum ^ {ni — 712} . Finally, by a rotation of coordinates: 



71 



r] = -^{xi- X2) , 



(119) 



Xl 



d 



d 



drj d 1 I' d _^ d 



dxi^2 dxi^2di dxi^2 drj -J2\d^ drj ^ 

also J2 becomes a semi-diffcrcncc of two isotropic harmonic oscillators: 



1 / & 



d 



1 



(120) 



and obviously the spectrum is the same. 

To analyze the decomposition of the infinite dimensional representation of 
SU{2) in irreducible ones we evaluate the Casimir J^. We start with 



J? = \ {P?Ql + P'iQl - [P1Q2P2Q1 + P2Q1P1Q2]} 

and from [Qa, Pb] = iSab, a, 6 = 1, 2, we have: 
1 



(121) 



7"^ 
"^1 



^ {PlQl + PiQl - [Pi (P2O2 + i)Qi + QiPi {Q2P2 - «)]} 

i {P2Q2 ^ p2Q2 _ ^p^p^Q^Q^ + Q1P1Q2P2) + Pi]} (122) 



while 



so that 



Then 



Ji = 4 {Pi Pi + Q\Ql + {P1P2Q1Q2 + QiQ2PiP2)\ 



jf + Jl = H{Pi,Qi)H{P2,* 



(123) 
(124) 



J' 



i [H{P,,Q,) - i7(P2, Q2)]' + H{P,, Qi)H{P2, Q2) ~ \ 
-^[H{P,,Q,)+H{P2,Q2)f ~\. 



(125) 



In the number representation H{Pa, Qa) Na + h, with A^i -I- iV2 = N, it is 



J' = l{Ni+N2 + lf 



4 2 V 2 



(126) 



This shows that the representation generated by Ji, J2, J3 contains all integer 
and semi-integer spins without degeneration, so it is a Schwinger representa- 
tion [G^. Finally, the basis ti)},— j < m < j, which diagonalizes {J^, J3} 
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is {\ni +n2,ni —712) = |7i+, 7i_)}, — n+ < n_ < ri+ and is obtained with op- 
erators N± = Ni zL N2 of the isotropic bidimensional harmonic osciUator, by 
rotating the old basis {|ni,n2) — \ni) \n2)} (see, e.g., 

The whole Hilbert space Ti decomposes in a direct sum of finite dimensional 
subspaces W of dimension 2j + 1, which are the carrier spaces of the spin 
tomography. In each Ti.^ , the {2j + 1)— dimensional spin tomography uses a 
set of fiducial operators, the eigenstates of the spin projection on z— axis J3, 
a generator of an irreducible unitary (2j + 1)— dimensional representation r of 
SU{2) : 

{Ao{j,m)},ne„ = {\jm){jm\},ne<T, cr = {-j, ,j} (127) 

The unitary transformations may be either the operators of the unitary 
group U{2j + 1) or simply those of the representation t of SU{2), depending 
only on the angle parameters which determine a point on the sphere : 

[/(<?), ffef/(2j + l); or C/(5) = r(g), g e SU{2). (128) 

In the last case the tomogram describing a spin state is a probability distribution 
depending on random discrete spin projection along the direction determined by 
the angle parameters. General reconstruction formulae for the spin tomography 
were obtained in [27]. 
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4 Quantum Mechanics in the tomographic pic- 
ture 



4.1 Superposition rule 

Having put the tomographic approach within a general setting of states S and 
observables O, we may now ask how one can add probabihties to describe the 
superposition of pure quantum states which allows for the description of in- 
terference phaenomena. For the case of two (orthonormalized) state vectors 
\ipi) and IV'2), given the probabilities (non-negative numbers) pi and p2, with 
Pi + P2 = Ij the superposition rule states that the linear combination 

IV') = VPr|V'i)+e*'"VP^IV'2> (129) 

is also a (pure) state vector. Here the factor e*"^ corresponds to the relative 
phase of the two state vectors. In the tomographic picture the formulation of 
the superposition principle should be discussed in terms of density states. In 
fact, since density states depend nonlinearly on state vectors: 

Pi = |^i)(V^i| , P2 = m{H, (130) 

a naive way to add them, that is pipi +P2P2, yields an operator corresponding 
to a mixed quantum state, different from a rank-one projector p = jV') (V'l 
corresponding to the pure state of Eq. (|129p . However, a way to add rank-one 
projectors to obtain a rank one projector has been proposed in |67j and is given 

by 

P = PlPl + vTlP2^======r , (131) 

VTr(piFoP2^o) 

where Pq is a fiducial rank-one projector corresponding to the arbitrary phase 
factor e^'^ of Eq. (|129p . The previous formula can be viewed as a purification 
of the mixed quantum state pipi + P2P2 ■ If Pi and p2 are not orthogonal, the 
formula yields a non-normalized p, so it requires a normalization factor (Trp) ^ . 
For the denominator to be different from zero we need to superpose only states 
which are not orthogonal to Pq. 

The composition rule for tomographic probabilities corresponding to the 
superposition of density states (Eq. ()13ip ) reads, using x = {X, fi, v) 



Tp{x) ^ piTp,{x) +p2Tp^{x) + ^==£^L=={Tp,p^p^{x) +Tp^P„p,{x)). (132) 

VTr(piFoP2^o) 

This equations contains tomograms such as Tp^p^^p^(x) which is possible to ex- 
press in terms of tomograms of pi and p2 , by means of the star-product formal- 
ism. 

In a star-product procedure, operators are replaced by their symbols, which 
are functions along the lines of section [321 

A^jAix). (133) 
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The associative product of operators is mapped onto an associative product of 
symbols 

AB ^ fAB{x) := fA{x) * !b{.S) . (134) 
This star-product is described through a nonlocal integral kernel as 

fAB{x) J fA{xi)K{xi,X2,x)fB{x2)dXidX2 ■ (135) 

Tomographic probabilities (tomograms) can be considered as symbols of opera- 
tors in a specific kind of star-product scheme. 

In [281 68 it is shown that the symplectic tomogram 7^ (X, /i, ly) is the symbol 
of an operator A in a specific instance of star-product, with star-product kernel 
depending on the continuous variables {X,ii, v) — x. It is given by I69j: 

K{Xi,^ii, 1^1,^2, Ai2, V2,X,^l. v) ^ -^5{^{vi + V2) - v{^i + /i2)) (136) 

vi^2 - i^2Mi + 2ai + 2X2 1 X 



X exp <; - 



Thus, the addition rule for tomographic probabilities, corresponding to the 
superposition of density states Eq. (jl3ip . reads, using again x = (X,fj,,v) 

Tp{x) = piTp^{x) + p2Tp^{x) (137) 

, Tp^ {x) * To (f ) * Tp,,_ {x) + Tp^ [x) *Tq{x)* Tp^ (x) 

"^^'^^ ^C{TpAS)*%{x)*Tp,{x)*%{x)) 

where, restoring [X, fi, v) : 

Tr(A) = C{Ta{X, fi, v)) = J Ta(X, fi, v)e'^ 5{^i)5{v)dXdpLdv . (138) 

We can describe the addition rule of Eq. (|137p in the following terms. Given 
the tomograms and Tp^ of the pure states pi, P2, and the tomogram Tq 
of a fiducial chosen pure state Pq, then the tomogram Tp is the tomogram of 
the pure state p — of Ea. (|131[) corresponding to the linear superposi- 

tion of Eq. (|129p . The formula Eq. (|137p guarantees that if 7^^ , 7^^ and Tq are 
tomographic probabilities of rank-one projectors, the result of the nonlinear 
summation Tp is again a normalized tomographic probability corresponding to 
a superposition state that can be realized in nature. 

The composition rule Eq. (jl37p has its partner in terms of Wigner functions 
of pure states. Let the Wigner function yVi(p, q) corresponding to a pure state 
\ipi) ("01 1 be composed with the Wigner function W2{p,q) corresponding to an- 
other pure state IV'2) (V'al and let the composition W{p,q) correspond to the 
superposition vector \tp) = y^lV'i) + ^^"^ \/P2\''p2) with the parameters used in 
Eq. p37p . Then one has [70], dropping the arguments {p,q), the result: 

W ^ PiWi + P2W2 + y/prp^ — 139 

y/C{Wi * Wo * W2 * Wo) 
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where, restoring the arguments {p,q), for any operator A and its Weyl symbol 
q) one has 

Tr{A) = <i)) = ^J^ <l)dpd<l ■ (140) 

The star product of two Weyl symbols is given by the formula 
Wi{p, q)*Wo{p, q) 
= j yVi{pi,qi)WoiP2,q2)K{pi,qi,p2,q2,p,q)dpidqidp2dq2 (141) 

with the Gronewold kernel 

K{pi, qi,P2, q2,P, 9) = ^ exp [2i{qpi - qip + qip2 - q2Pi + q2P - P2q)] ■ 

One can easily sec that the exponent in the r.h.s. of the above formula may be 
written as 4iS', where 5 is the area of a triangle in the phase space with vertices 
in the points (91,^1), (92,^2), {q,p)- 



4.2 Uncertainty relations 

Having considered the superposition rule within the tomographic scheme we 
consider now the formulation of uncertainty relations. We present below the 
derivation of uncertainty relations in a general framework. Given an operator 
A, then it is obvious that 



{A^A) = Tr (pA^A) > 



(142) 



for any density state p which is a non-negative hermitian operator with unity 
trace. Let A be, for systems with one degree of freedom, a linear combination 
of operators {Bk} to be suitably chosen. 



A 



Then 



(143) 



(144) 



that can be rewritten, dropping the label p and using commutator and anti- 
commutator, as 



> 0. 



(145) 



Since c^, Ck are arbitrary complex numbers, the positivity of the quadratic form 
means the positivity of the matrix 



iS) 



hk 



Bl,Bk 



Bl,Bk 



(146) 
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The Sylvester criterion provides the necessary and sufficient condition for the 
positivity of the matrix in terms of its principal minors 



Mi{S) > , A'hiS) > , 



det 5 > . 



(147) 



This general scheme allows to obtain uncertainty relations. 

Thus, for position Q and momentum P we choose the error operators as B's 



Bi = Q - (Q> - AQ , B2 = P-{P)= AP 



So, we get 



S 



(^QP CFPP 



{[Q,Pl 
{[Q,P]) 



where we have introduced the variances 

\2 



<^QQ = (0') - {QY , ^pp = {p^) - {py 



and the covariance 



(148) 



(149) 



(150) 



(151) 



We then obtain from the positivity of the Mi principal minors the obvious 
conditions 

<yQQ>0 , CTPP > , (152) 
and from the positivity of the determinant 



det S - aggapp - [ a^p + - {- [Q, P])' ] > 



the Schrodinger-Robertson uncertainty relation follows 

{AQ^){AP^)>l{[AQ,AP]l^ , ^, 



(153) 



(154) 



this reduces to the weaker Heisenberg uncertainty relation for uncorrelated 
states * 

{{P') - (Pf) {{Q') - (Qf) ^ \ , {n = i). (155) 

In view of Eq. ([ST|) . we can write Eq. p55p as: 



J X^T{X,0,l)dX - J XT{X,0,l)dX 
J X^T{X,l,0)dX - J XT{X,l,0)dX 



(156) 



1 

> - 

- 4 
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because T(X, 0,1) is the momentum probability density and T{X,1,0) is the 
position probabihty density. The uncertainty relation in covariant form reads 



X^T{X, cos9,sm9)dX - 



X^T{X,-sm9,cose)dX 



J XT{X, cos6',sin6')dX 
J XT{X,-ain0,cose)dX 



(157) 
1 



> 
- 4 



The obtained relation takes into account all the values of optical tomogram for 
all angles 9. It allows to use the values of experimentally measured optical to- 
mogram to check the uncertainty relation of position and momentum by using 
a homodyne detector. This formula, that can be written in the short form 
> 1/4, provides a constraint on admissible quantum mechanical tomo- 
graphic probabilities due to Heisenberg uncertainty relations. The experimental 
check of the Heisenberg uncertainty relations can be done using the positivity 
condition of the function (uncertainty function) 

$(6») := F{9) - ^ > . (158) 

Such a checking can be done [7T] (see also the discussion in [72]) using experi- 
mental data of [n nn m [m . 

The tomographic expression for Schrodinger- Robertson relation requires the 
star product to represent the mean value of the anti-commutator. 

For the case of angular momentum we choose selfadjoint B's operators as 

Bk^AJk , A: = 1,2,3 (159) 

and get, recalling that [Jh, Jk] = i^hkiJi ■ 



<Jj,j, <J.hJ. +o -(-^3) (Ji) (160) 




where again variances and covariances appear. 

The three first principal minors give the obvious conditions of positivity of 
all the variances 

CTjiJi > , cr,/2j2 > , fjjgjg > 0. (161) 

The first of the three second principal minors yields a quadratic relation with 
indices 123 

^J^J^^J.J. - (^,A,7, + \ {Jzf^ > (162) 

while the others give quadratic relations obtained from the first by circular 
permutation of the indices 231 and 312. 

Finally, the positivity of det S yields a cubic relation connecting the vari- 
ances and covariances of all the components of the angular momentum, or spin 
projections. The tomographic expression for a star-product kernel providing 
such a cubic relation exists, even though cumbersome [75] . 
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4.3 Classical and quantum distributions: examples 



Some observations are now in order to discuss the conditions the uncertainty 
relations give for considering a 'tomographic' normalized function the quantum 
tomogram of a density state (the necessary and sufficient conditions for a func- 
tion on phase space to be a Wigner function are discussed, e.g, in [7B]). As we 
have observed previously, an exact condition stems out from the positivity of 
density states. Thus, given a 'tomographic' probability T{x) (either for spin 
variables or for continuous ones as position and momentum), the appropriate 
reconstruction formula can be used to get an operator gr- Then, we can check 
wether its eigenvalues are non-negative. If they are, gr is a density state, pos- 
sibly non-normalized, and T is a quantum tomogram. However, this condition 
is rather formal and other, possibly more operative, conditions may be given. 
The 'tomogram' T is assumed to satisfy the specific properties which hold for 
the tomographic probabilities in the given tomographic scheme, such as the 
homogeneity condition T{XX, X^, Xv) = T(X, \X\ in the symplectic case, 

for instance. Altogether these properties guarantee the hermiticity of gx and 
the normalization Tr^gq-) — 1, but they are not sufficient for the nonnegativity 
of gT- For that, there are necessary, in general not sufficient, conditions like 
uncertainty relations. In fact, the nonnegativity of the operator gr guarantees 
the fulfilling of all available uncertainty relations because they stem from the 
equation (AM)^^ = Ti{A'^ Agq-) /TT{gr) > . Of course a complete (usu- 
ally infinite) set of uncertainty inequalities is sufficient for the nonnegativity of 
gT- However, such a completeness is a tautology, as in a sense it stands for 
'positivity'. 

So, if the corresponding tomographic expressions of such inequalities are 
satisfied by T, we have the necessary, or sufficient, conditions for T to be a 
genuine quantum tomogram. 

This method requires the use of the dual tomographic map [77^ for the 
operator symbols f^ix) :— {A\G{x)), as for an observable A one has |78j : 



{A),^ = TiiAgr) - {A\gT) = / {A\G{x))(P{x)\gT)dx ^ / T{x)fi{x)dx 



Thus, the Schrodinger- Robertson uncertainty relation can be written for the 
symplectic case as: 




(163) 
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where {Q) , (P) are assumed to be zero, while the dual symbols are given by: 



{X, fi, v) = ^Tr(Q2 exp [i(X - fiQ ~ ,,P)] , (164) 

,f$. (X, ly) = -^Tr(p2 exp [i{X - fiQ ~ lyP)] , 
Zn 

/|[Q,P]+ = ^Tr(i[Q, P]+ exp [i{X - /iQ - z.P)]), 

and are known generalized functions [TSl [HO] ■ The Schrodinger-Robertson re- 
lation, written in terms of these generalized functions integrated with the T 
function, provides a necessary condition for T to be a quantum tomogram. We 
illustrate this analysis with the help of some examples. 
Example 1 Let us consider the probability distribution 



and check wether it fulfills or not the uncertainty relations. 
To this aim one has to evaluate dual symbols like 

f^.{X,f,,,y) = i-TV(02e^(^-A<Q--P)). (165) 
By using the Weyl symbols 

A^WA{p,q) (166) 

the trace is readily evaluated as 

Tr(AB) = / ^^^^^^dpdq. (167) 



We have 



so that 



Wq2 = q\ (168) 



f^X,,,.) ^ ^Tr(QV(--^--)) = ^/l^^- Idvdq 

^ dq = -e'^ 6{iy)S" (fi) (169) 

27r 

and analogously 

fUX,f^,^) - -e'''S"iiy)Siti), (170) 
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Thus we obtain 



2[AQ,APl 



v/27r(/i2 + ,.2 



-.e'^5{v)"5(pL)dXdiidv = - 
I 2 

-.e'^5'{v)5'{^i)dXd^idv = 



so the Schrodinger-Robertson uncertainty relation 

(AQ2)(Ap2)>l([AQ,AP]+)Vi 



(172) 



is satisfied. In this case it is easy to reconstruct the quantum state qt^ it is the 
ground state of a harmonic oscillator. This can be done or directly, by the use 
of Eq. or in two steps, first by a Radon anti-transform of the tomogram 

yielding a Wigner function on phase space, from which the quantum state is 
readily obtained. 

Example 2 An example of a non-quantum, classical tomographic probability 
distribution is provided by the following normalized distribution 



1 



+ 



where I{Q,ijl + v) is the interval with extrema and fj, + i>. As it is homogeneous, 
it has the form of a symplectic tomographic distribution. We have 



(Ag^) 



1 



X/CCm+i') {X)e'^ S{i')5" {fi)dXdfidv 



+ 



■Xiio,MXyS"{fi)dXdfi 



/■ 2 sin ^ 


■ M" 


/ exp 




/ 


L 2J 



S"{n)dfi 



Analogously 



\ ' 3 



(173) 

_ 4 
~ 3 

(174) 



While 



[Ag,AP]_ 



2siniii^ 



l/i + v\ 



XHo.,.+.){X)e'''5'{,y)S'ip)dXdf,d,y 



exp 



S' {v)6' ifi)dtid:^ = - 
o 



so the uncertainty relations are violated: 



4 4 
3'3 



(175) 



(176) 
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To obtain the classical distribution, we have to evaluate the Radon anti- 
transform of the tomographic distribution : 



1 



1 



(27r) 

(2n? I 



2sin(i^) 



Xi(o,iJ.+v){X) exp [i{X - iiq- vp)] dXcifidu 



(27r)" J IJ' + J^ 
By introducing the new variables 



exp 



- Iiq-up 



diidv 



we then get 



4 = 



sin^ 



^ — v 



11 = £,+iT,i^ = £,-r] 



(177) 



(178) 



(27r)- 



exp - ^{q + p)- rjiq - p)}] d^drj 



(179) 



sin^ 



exp[i{^-^{q+p)}]d^ 



I 

= X[-i,i](i- iQ + p))HQ-p) = X[o,i]{q)HQ - p) 

Example 3 Another non-quantum, classical example is given by the positive 
homogeneous normalized distribution 



1 



2|Ai + ^| 



Remembering that 



1 



2\n + i' 



exp 



exp 



\X\ 



\X\ 



iM + i^l 



l/^ + i^l 



e'^dX = 



l + \n + iy[ 



we get 



{AQ^) = {AP') = 2;(-[AQ,AP] 



(180) 



(181) 



so that again the uncertainty relations are violated, and the tomographic dis- 
tribution is not quantum. With the change of variable 



in the Radon anti-transform, it corresponds to the classical distribution: 



(182) 



1 



(27r)^y 1 + lM + i^l' 

^^/iTi^'^P 

^^(^(a-p))fexp 



■ exp \—i(jiq -|- vp)\ d\xdv 
-i\{(l^P)-^\{(l-p) 



d^dr) 



\q + p\ 



1 



%-p)exp[-|g|] (183) 
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Example 4 The analogous approach to spin tomography is based on the neces- 
sary conditions of Eq.s (|1611I162"1) and positivity of the determinant of the matrix 
S of Eq. (|160p . All the elements of S can be written in terms of spin tomograms 
and dual spin tomographic symbols of the observables Jk, JhJk{k,h = 1,2,3), 
which can be obtained using the Gram-Schmidt operators for the spin tomo- 
graphic star-product scheme given in [75]. For instance, consider in the qubit 
tomographic case the vector component 



9 9 1 1 

Ta,l3 in)^acos'^-+ (3 sin^- = - {a + (3) + - {a - P) cos9, 



the other component being 1 — Tajj i^) , corresponding to the operator 

Pa.,P = 



a 
/3 



The elements 



1 



of the S matrix, in view of 



JhJk + JkJh 



(184) 

(185) 
(186) 

(187) 



read as 

dQ {5Hkh\Gn)T^,p m 



with 



while 



dn{A\Gn)To,,p{n) I dn{Jk\Gn)T^,pm 

(188) 



1 



dfl (ShkhlGn) Ta.,p {^) = -{a + P) 5hk, 



{•h\Gn) 
{J2\Gn) 
{■h\Gn) 



2ti 
3 

3 

47r 



■ cos sm c 



{J2) 









008 6* =4> (J3) = - (a - /3) 



so that eventually 
5^1 



[a + (3) i{a- 13) 
-I (a -(3) (a + (3) 








[a + 13) -{a- (3Y 



(189) 
(190) 
(191) 

(192) 



Now S is non-negative iff < a, /? < 1. Then, < 7^,/3 (fJ) < 1 is the component 
of a probability vector corresponding to the (non- normalized) density state Pa,j3- 
When the uncertainty relations are violated, Ta^js (fi) is readily recognizable not 
to be a probability vector. 
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4.4 Equations of motion 



For the tomographic description of quantum mechanics to be equivalent to other 
existing ones we have to consider the formulation of dynamical evolution, along 
with the integration of the equations of motion usually done by solving an 
eigenvalue problem. In the coming two subsections we are going to address 
these aspects. 



4.4.1 Time evolution 

The evolution of a state probability distribution in classical hamiltonian me- 
chanics with potential U (g) is given by the Liouville equation 

rj r\ fi 

g-^f(l'P^t) +Pg^f{'l,P,t) - Q-U{q)—f{q,p,t) = (193) 

We study now evolution in quantum mechanics using Wigner functions and 
tomograms. We recall the definition of Wigner function of a density matrix 

p{q, q') ■■ 

P(9 + 2 ' 9 - 2 ^ exp(-ipa;)rfa; . (194) 

The equations of motion for p, derived from the Schrodinger equations, give rise 
to equations of motion for yV(p, q). 

The evolution equation for the density matrix p{q, q',t) is obtained from the 
Schrodinger equation for the wave function. In fact, for the Hamiltonian 

H{P,Q) = ^P^ + U{Qy,{m = l) (195) 
one has in the position representation: 

and 

-ilr{<l',t)= (^-l-^ + Uiq')^r{q',t). (197) 

So, upon subtracting the above equations after multiplication by 'ip*{q',t) and 
ip{q, t) respectively, we get the von Neumann evolution equation for the density 
matrix p{q,q',t) = ip{q,t)^p* {q' ,t) of the pure state tp : 



i^PiQ,<l',t) 



1 ff^ 1 ff^ 



p{q,q',t). (198) 



Equations of motion for Wigner function yV{p, q) can be derived by using the 
following association among derivatives of the density matrix and the Wigner 
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function: 



l-^W{q,p) 
qW{q,p) 



9 d , , 



\ {q + q)p{q, q) 
\ f d d 



(199) 



and 



qp{q,q') 
q'p{q, q') 



1 d 

7:^-w]yv{q,p) 



ld_ 



(200) 



Then, using the previous correspondence rules, we finaUy obtain the evolution 
equation for the Wigner function yV(g,p, t), which reads 



i^W{q,p,t) 



ld_ 
28^ 



ip 



ld_ 
2d^ 



+ 



U q 



2 dp 



-U[q 



ip 

i_d_ 
2 dp 



Wiq,p,t) 
W{q,p,t) . (201) 



The Wigner function is real. The above equation can be written as 



d_ 

dt 



+ ) W(g,p,t) 



U[q 



i_d_ 
2dp 



U q 



i_d_ 
2dp 



W{q,p,t) 



(202) 

These equations of motion may be called the Moyal evolution equations. The 
expansion of Eq. (|202p up to first order in (^^d/dp) yields the classical Liouville 
equation p93p 

A similar correspondence table allows us to derive the evolution equations 
for tomograms out of the evolution equation for the Wigner function. To get 
the evolution equation for the tomogram Tp{X, fiTV^t) the correspondence table 
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which is obtained from the Radon transform formulae is 



d d 
— ^ u—T{X,i^,u) 



d f d 



(203) 



Thus, from the Moyal equation we obtain 

d d 



T{X,iJi,u) 



dt ''dv 



T{X,ii,v,t) = -i 



U 



dn \dx) 



— V 

2 dX 



c.c. 



T{X,ii,v,t). 
(204) 



The "Planck constant expansion" of both Moyal and tomographic evolution 
equations is given by the potential energy expansion in powers of (— ^5/9p) for 
Moyal and [—^vd/dX^ for tomographic equation, respectively. For example, 
the Moyal equation for an harmonic oscillator is 



|+p|)>V(,,p,i) 



while the tomographic one reads 



|-M|;)T(X,M,.,i) 



i d 



c.c. 



yV{q,p,t) 



d f d 



djji \dX 



i d 

+rdx 



c.c. 



-iy—T{X,iJ,,u, t) 



(205) 

T{X,ii,y,t) 
(206) 



4.4.2 Eigenvalue problem 

We may obtain the Moyal and the tomographic form of the eigenvalue equation 
if we start from the von Neumann stationary equation for the density matrix. 
In the position representation, we have 



Ep{q,q') = 



1 92 



1 92 



' 2 5g2 2 dq' 



■ U{q) + U{(i) 



p{q,q'). (207) 
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so that with the previous correspondence rules the Moyal equation for the 
Wigner function is 



EWiq,p) - -- 



1 d 



1 d 

2 9^-*^ 



rr I i d \ ^ / id 



W{q,p) 

W{q,p). (208) 



and the corresponding tomographic form of the equation: 



\ d_ ■d_ f_d_ 



+ c.c. 



U 



a / a 



dfi yax 



i a \ 

raxr'-'- 



T{X,tiAW 



for the energy levels E. In the usual pictures, the energy spectrum is obtained 
by solving an eigenvalue problem of the Hamiltonian operator in the carrier 
Hilbert space. Thus, for the n— th level, Eq. (I203|) is to give the result 



Tn{X,n,v) 



-Hi 



X 



(210) 



However, neither the Wigner function W{q,p) nor the tomographic proba- 
bility densities T(X, /i, v) are vectors of a linear space. The eigenvalue problem 
in the Hilbert space of states is mapped onto the problem of solving integro- 
differential equations either of Moyal (von Neumann) form or equations for 
probability distributions. It is remarkable that in time evolution ruled by the 
tomographic equation, the function T{X, /i, v, t) at any time t is positive and 
normalized if the initial tomogram is. This is obvious due to the correspon- 
dence chain rules substituted in the initial unitary evolution of von Neumann 
equation. Nevertheless it would be difficult to see it at a first glance directly in 
the tomographic evolution equation, without going back to the von Neumann 
equation. 



4.5 Composite systems: separability and entanglement 

To define the notion of separability of quantum state of a composite system let 
us discuss the example of two qubits. Each of the qubits has in its own Hilbert 
space of states the standard basis {|to)} where m — ±1/2 is the spin projection 
on the z-axis, that is 




(211) 
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The spin operator S = {S^, Sy,Sz) is determined by the PauU matrices {a^, (Ty,az) 
as 

(212) 



S=-a {h=l) 



where 



1 

1 



-i 

1 



1 
-1 



Thus the basis {|m)} satisfies the eigenvalue problem 

Czlm) = m\m) . 



(213) 



(214) 



For two qubits the Hilbert space of states is four-dimensional and it is ob- 
tained as tensor product of the two-dimensional Hilbert spaces of each qubit, so 
that we have four basis vectors 



|?Tiim2) = |TOi)|m2) 
The previous equation means 



mi,m2 = ±^ 



(215) 



1 1 

2^2 



1 1 

-2 + 2 



1 1 

2 ~ 2 






V J 

1 



V / 



1 

V / 




V 1 J 



(216) 



(217) 



(218) 



(219) 



Any pure state of one qubit \tjj) can be expressed as superposition of the basis 
states 



tp2 



The density matrix for one qubit in case of a pure state reads 



(220) 



(221) 
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Pll Pl2 
P2I P22 



The density matrix for a mixed state of one qubit 

P 

is such that 
and 



P = P^ 



Tip = 1 



(222) 

(223) 
(224) 



p > 4=> Pll > ; P22 > ; P11P22 - P12P21 > . 
It can can be presented in the form 

p = Ai|Vi)(V'i|+A2|^2)(V'2| , (225) 

where Ai,A2 > are its (non-negative) eigenvalues and |V'i):|'02) the corre- 
sponding eigenvectors. 

The simply separable state of two qubits is defined as the factorized state 



IV) = \'Pl)\f2) = 



mi 

^12 



V21 

f22 



<y5ll'^22 
<y5l2'^21 
V <y'l2'P22 ) 



(226) 



The density matrix of this state is the tensor product of the density states of 
each qubit: 

Pip = W) (<^l = Pvi ® Pv2 ■ (227) 
A separable state of two qubits is defined as a convex sum of the above factorized 

StcLtCS 

p(l,2)=^^AfepW®pW . (228) 
In Eq. (|228p the states , plp2 satisfy the condition of purity 

{pI%) {pI'L) - pI% ■ (229) 

However, this condition can be violated, as the general definition of separable 
two qubit state is: a state is separable iff its density matrix can be given the 
form 

~ (230) 



P(1'2)=V A 



(fe) 



kPi 



where X]fc Afc = 1 with Afe > 0, while p'"^\ p^2^ are density matrices of each qubit 
respectively. 

The above separability condition is written for density states. A condition 
which holds for the corresponding tomograms may be obtained by calculating 
the tomograms of the two sides of the above equation. The tomogram of the 
state p(l, 2) is the joint probability function: 



•P2 



(fc) 



Tpii,2){'mi,m2,U) = (mim2 \Up{l,2)U^\ mim2) 



(231) 
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As the tomogram of a simple separable state factorizes in the product of inde- 
pendent probabilities: 



TOim2 
mi 



Ui (g) U2 {pi ® P2) ul ® ul 



UipiUi 



nil) (1712 



U2P2UI 



77117712 
7712 



(232) 



we eventually get from Eq. (j230p the tomographic separability condition in the 
form: 



r,(i,2)(7ni, m2, Ul ® U2) = J2 AferW("^i, c/i)rW(™2, c/2) 



(233) 



The state is called entangled when the 4 x 4-density matrix p(l,2) cannot 
be presented in the form (|230p . or its tomogram 7^(1 2) in the form of Eq. (|233p . 

Now we reformulate the introduced notions and definitions using the tomo- 
graphic representations of the qubit states. This representation is characterized 
by a map from density matrices onto the family of the probability distributions 
which is invertible. We racall that for one qubit one has 



Pii P12 
P21 P22 



Tp(+i 71) 

M-2'^) 



(234) 



where the unit vector 71 = (sin 6 cos (p, sin 9 sin ip, cos 9) determines a point on the 
sphere S'^. As it was observed in section IXSl the tomogram here is presented in 
the form of a probability vector, whose components Tp{m, n) = Tr(U^ \m) (m| Up) 
are the diagonal elements of the unitarily rotated density matrix: 



M+^,n) = {UpU^)n ; Tpi-^,n)^{UpU^) 



22 



where 



U 





Ul2 \ 




V ^t21 


U22 ) 





COS ^e'('>+v)/2 



sinfe^('>-^)/2 
sin|e-^('^-^)/2 cos f e-*(^+^)/2 



(235) 



(236) 



is a unitary matrix of the group SU{2) parametrized by the usual Euler angles 
9, tjj, ip. The matrix U transforms the spinors according to the rotation labeled 
by the Euler angles. 

For the qubit the tomogram can be considered as a function on the group 
5/7(2) (in fact on the homogeneous space SU{2)/U{1)). The formulae can be 
presented in the form 



\un\ pii 

\2 , 



|Ml2pP22 + (ui2Un/021 + C.c) 



\U2\\^P\\ + |m22|^P22 + (M21U22/O12 + C.c). 

In terms of Euler angles the tomographic probability reads 



(237) 



cos^ -pii + sin^ -P22 + cos - sin -(e "^p2i + e"^pi2) 



Tp(--,C/) = sin2-pn+cos2-p22-cos-sin-(e*'^pi2+e-"^P2i)(238) 
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One can regard the qubit tomogram in the following manner. The density 
matrix ip can be considered as a point in the Lie algebra u{2) of the group U{2); 
then the formula (|235p defines an orbit of the unitary group in its Lie algebra. 

Assume pi2 = P21 = 0: the initial point of the orbit, i.e. U = /, is determined 
by the two non-negative numbers pn and P22 which satisfy the simplex condition 
Pii + P22 — 1, that is a segment with extremes in (1,0), (0, 1) in R^. Thus we 
may choose in the Lie algebra of the group U{2) a subset, which is a simplex 
in an afSne subspace modeled on the Cartan subalgebra of SU{2). The formula 
for the tomogram (|237p determines the orbit of the group SU{2) in the simplex. 
Further details of this relation between qudit tomograms and points of a simplex 
are discussed in [ST]. It is interesting to note that, under the action of the group, 
the orbit starting from any point of the simplex does not go out of the simplex. 
Thus one can define the tomographic probability as a map of the pairs {p, U) 
onto the points of the simplex Tp{U). Here p = {pii,p22)'^ is a column vector 
and U is an element of the SU (2) group. As we remarked in section 13.91 we 
could also use the full unitary group in Eq. (|237p . 

The general case pi2 = P21 7^ is analogous due to the possibility of diago- 
nalizing: 

Ul,U„ ^ («■ I) . (239) 

so the generic state p is again mapped onto the orbit of the group in the given 
simplex, but the initial point of the unitary transformation is shifted by the 
diagonalizing unitary matrix Uq, i.e. now the element UUq moves the simplex 
point (pii,p22)- 

We stress that the initial points of the orbits belonging to the simplex can be 
considered as elements of a Cartan subalgebra of the unitary group belonging 
to non-negative weights (the Weyl chamber from which all other weights can 
be obtained by discrete reflections (Weyl group)). Thus the formula for the 
tomogram (I237p can be written as: 

Tp{+^,U) = |(f/C/o)ii|'pii + |(C/C/o)i2pP22 

Tpi-l^U) = |(t/C/o)2i|'pii + |(C/C/o)22pP22 , (240) 

in the form of a bi-stochastic map acting on the simplex where the shift ma- 
trix Uo and simplex point coordinates (^11,^22) are connected with eigenvalues 
and eigenvectors of the density matrix p by Eq. (j239|l (for further details, see 
[51]). The equation (|240p can be reinterpreted in the following way. The ortho- 
stochastic 2 X 2— matrices 



M{U) = 




(241) 



belong to the semigroup of bistochastic 2 x 2— matrices. 

We recall that, in the 71— dimensional case, if e G M" denotes the column 
vector with all components +1 and its transpose, an n x n-matrix M is 
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called (column) stochastic iS e^M = and bistochastic iff 

Me = e and e^M = . 



(242) 



All such M's are non- negative matrices. For invertible matrices M, if we consider 
not only bistochastic matrices but also their inverse, which need not be non- 
negative any more, we get an open dense subset of GL{n — 1,M), the general 
linear group of {n—l) x (n— l)-matrices ^T]. Then formula (|240p . when U varies 
on all the elements of the unitary group, provides an orbit of the semigroup of 
matrices M(U) in the simplex. Thus the tomogram of the quantum qubit state 
is orbit of the semigroup in the simplex, which is a subset of the first positive 
chamber of the Cartan subalgebra of SU (2) group. The statement is correct 
for unitary tomogram of any qudit state. Since a tomogram in the tomographic 
representation is identified with a quantum state, for a qudit one can say that 
the quantum state is the orbit of the semigroup of bistochastic maps M{U), 
parametrized by the pairs {([/, 7^(C/))} which are the graph of the tomogram, 
in the simplex which is a subset of a chosen Cartan subalgebra of SU {n) group. 
Due to this picture we write the qudit analog of Eq. (|240p in matrix form as 

fpiU) = MiUUo)7. (243) 

Here p = (pii, P22, • ■ ■ , Pnn)'^ is the probability column vector consisting of the 
eigenvalues of p or, equivalently, a point in the simplex in the given subalgebra. 
The ortho-stochastic matrix M{UUo) has matrix elements 

{M{UU,))^^ = mUoW- (244) 

The columns of the matrix Uq are normalized eigenvectors of p, and a 'gauge' has 
been chosen by fixing the phase factors of Uq and an ordering of both the com- 
ponents of p and the columns of Uq so that UqpUq = diag [pii, P22, • ■ ■ , Pnn] ■ 

The component of the vector Tp(U) are tomographic probabilities. They are 
defining a spin tomogram if one takes as n x n-matrix U e SU{n) the matrix 
of a (2j + 1 = n) —dimensional irreducible representation of the group SU{2). 

For two qubits the condition determining a separable state of a composite 
(bipartite) system can be rewritten in the form 

fp(i,2) (f/i ^U2)^Y. ^'^ {pM'^^I' ® \U2Ui^^ pf ' ® pf ) (245) 

k 

or 

|(f/i <E> U2)Uo{l,2)\^p{l,2) = ^ Afc {\UM^^\^ ® \U2Uio''\') pf ^ ® ^ , (246) 

k 

where the notation \A\'^ means (|Ap)js — jAj^p. 

Thus, given the (non-negative) eigenvalues and eigenvectors of the density 
matrix of two qubit systems p(l,2), which are the components of the vector 



48 



p(l,2) and the corresponding columns of the unitary matrix Uo{l,2) respec- 
tively, the state is separable if the vector in the left hand side of Ea. (|246p is 
a convex sum of vectors, which may be written as tensor products of the to- 
mographic probability vectors of each qubit, i.e., with eigenvalues given by the 
components of pi , p2 and eigenvectors given by the columns of the unitary ma- 
trices u[q^ and ■ 

In other words, the state is separable if the semigroup orbit determined 
by semigroup \{Ui (g) [/2)t^o(l, 2)^ acting on the simplex can be presented as 
a convex set of orbits of semigroups determined by semigroups |J7iC/}q''P and 

1^2^20^ P acting on their own simplex. 

Now we illustrate the previous theory by discussing some examples. 
Example 1 Let us consider the simplest example of simply separable state of 
two qubits: |tt) = |T)i 1 1)2 ■ The density matrix of this state is diagonal: 



P(l,2) 



/I 0\ 





yo 0/ 



(247) 



The eigenvectors of this matrix may be chosen as 






















1 


















;wo3 = 


1 



















vJ 



(248) 



and the eigenvalues are (1, 0, 0, 0),.that is: 

^(1,2) = 








(249) 



The unitary matrix Uq constructed from the eigenvectors is, of course, the iden- 
tity matrix: 

/I 0\ 

c^o(i,2)- 0010 
\o 1/ 

The tomogram Eq. ()246p of the state reads 



because 



fp^iMUi U2) = |c/i®f/2|V(i,2) 

(Ui (g U2)Uo = {Ui ® U2)h = Ui®U2 



(250) 



(251) 
(252) 
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and thus Eq. (l25ip yields 

^p(l,2)(t^l'»t^2) 



/cos^ 


2 


cos^ 


%\ 


cos^ 


ei 


sin^ 




sin^ 




cos^ 


k 


\sin^ 


I 

2 


sin^ 


2 / 



- |C/iC/io|Vi'^|t^2C/2o|'p2 (253) 



Here: J/n 



P 1 P2 



The formula (|245p contains only one term, with = 1 times the tensor 
product of the 2— vectors 



|C/iC/io|Vi = 



cos^ ^ 



; |«72i72o|V2 = 



cos^ f 
sin^^ 



(254) 



Example 2 Another example of simply separable state of two qubits is the state 
with both spin-projections oriented along the x-axis. The state vector reads 



V2 VI 



1/ vo 



1 /I 1 



^/2 VI 



1 VO 



1 
1 

VJ 



(255) 



Thus one has the density matrix of the two qubits state, which is a pure simple 
separable state: 

/I 1 1 1\ 



P(l,2) = - 



1111 
1111 

Vi 1 1 1/ 

The eigenvalues of this matrix yield the vector 



(256) 



P(l,2) 

The corresponding eigenvectors of the matrix p(l, 2) may be fixed as 








(257) 







('] 










1 




-1 




1 




-1 


1 


;«02-2 


1 




-1 


; "04 = 2 


-1 


vv 















The above vectors are the columns of the unitary matrix Uq: 



(258) 



/I 1 



1 



1\ 



I- 11-1 

II- 1-1 

VI -1-1 1/ 



^2(1 -1 



^2(1 -1 



(259) 



50 



The matrix {Ui (E) U2)Ua is a product: 



Here the matrices Ui and U2 have the form of Eq. (I236P 
The corresponding orthostochastic matrix is 

|(f/l®C/2)C/oP 



due to the obvious property 



\a^b\^ = \a\^(g,\b\\ 



(260) 



(261) 



(262) 



^From this orthostochastic matrix and the vector (|257p we get the tomogram 
probabihty vector as: 



Tp(i,2)(t/l®C/2) = |((7i®[/2)[/o|V(l,2) 

|2 



(263) 



/ 



VI 



^e^-^i + sin ^fe-^v'i I 



Icos^e*"^! + 
|-sin|Le'¥'i 
-sin^e^'^i 4 



2 

COS ■ 



COS ^e'-f^ H 



■ sin 



«V2 



|cos%e''^2 



sin^e''^^ 



sin^e-^'^^l 
f cos%e-*'^2|V 



The tomogram has the form of tensor product, corresponding to only one term 
with Afc = 1 in Eq. ([M5l) : 




(264) 



|cos^e"^i +sin|Le~*'' 
sin ^e^'^i + cos ^e-'^^^ 




cos ^e'^P^ 



sin %e^'^2 + cos ^e'"^^ I 



Finally, taking a convex sum of the density matrices (|247p and (|256p : 



p(l,2) = cos^S 



(I 














1 


1 


1\ 
















1 


1 


1 


1 














+ - sin^ 5 
4 


1 


1 


1 


1 


Vo 








0^ 




^1 


1 


1 


1/ 



(265) 



we get a separable state whose tomogram, by construction, has the form of the 
convex sum of Eq. (|245p . with Ai = cos^ 5, A2 ~ sin^ 5 : 

/cos2 |l cos2 ^\ 



I 



cos"^ ^' s\Y? Y 

sin^ ^ cos2 ^ 
Vsin2|Lsin^f / 



sin^ 5 



sm 



Vl 



|cos ^e^"^! + sin ^e-'f^\ ■ |cos f-e*'^^ + 
cos ^e*"^! + sin ^^e^^'^i | • | - sin ^e^^^ + cos %e 

- sin ^e*'^! + cos ^e'^-'^^f ■ |cos %e*'^= + sin %e 

- sin ^e'f^ + cos ^e^^'^i I • I - sin ^e"^^ + cos ^e'^'^^ I V 



(266) 
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ExEimple 3 Let us now consider an example of an entangled two qubit state 

(267) 



1 

71 



The density matrix reads 



IV') (V-l = 



/o 








1 1 
1 1 








0/ 



the four eigenvalues yield the probability vector 



P = 







The corresponding eigenvectors may be chosen as 



/1\ 





,U02 



1 

71 



1 

1 

voy 



/o\ 

1 

-1 

voy 






VV 



and the unitary matrix Uq diagonalizing the density matrix is 



\U01,U02,U03,U04,\ 



10 

^ ^ 

^ -4= 



Vo 



V2 





V2 







ly 



The second column of the matrix UUq has the form: 



1 
72 



I' Ui2 + Ui3 ^ 
U22 + U23 
W32 + W33 
\ U42 + U43 / 



(268) 



(269) 



(270) 



(271) 



(272) 



Then the second column of the corresponding orthostochastic matrix {M{UUo))-^ = 
I (£/[/o)js P is the vector 



P = 



\Ul2 
\U22 ■ 
\U32 ■ 



■Ui3\ 
■■"23 I 
■"33 I 



(273) 



\ \U42+U43\ I 
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Thus, the tomogram of the entangled two qubit state (|267p is the probabihty 
vector 

/ \U12 +U13P \ 



\U22 + U23\ 
\U32 + 
\ |U42 J 



For the subgroup U = Ui ® U2, where 



\ ai21 022 

the tomogram of the state reads: 

fpiU - C/i ® C/2) = ^ 



U2 



bu bi2 
621 622 



/ |aii5i2 + 012611^ \ 
\a11b22 + 012^21^ 



|a2i^i2 + a22fciiP 
V 1021^22 + 022^21 P J 

By using the Euler angles as parameters: 



U2 



an ai2 

0-121 0-22 

b\i bi2 

6121 ^22 



sm -jc" 



cos^e^('''=+'^^)/2 



sin^e*('^2-^2)/2 



-i(l/'2+V2)/2 



(274) 



(275) 



(276) 



(277) 



the components of the above tomogram are the following explicit functions of 
(ni,n2) e 52 (g) S"^ : 

1, 



^ , 1 1 ^ ^ , 
^p(+o> + o>"i:"2) = 



2 Oi . 2 C2 

„,cos — sm — 
2' 2 2 



• 2 "^l 2 "^2 

sm — cos — 
2 2 



(71 . (71 02 . 

+ COS — sm — cos — sm „ 
2 2 2 2 



(278) 



1 



:,ni,n2) 



i 2 (71 2 '72 

- cos — cos — 

2^ 2 2 



• 2 "71 . 2 ^72 

sm — sm — 



(71 . (71 

COS — sm — 
2 2 



cos — sin— (e'^'^i-'^^) 



^)] 



1 



2' ■ 2'"^'"^) 
-,--,ni,r?2) = Tp( 



1 1 



:,"i,"2) 



2' ■ 2' 

The problem of separability amounts to the existence of a decomposition of the 
above probability vector as a convex sum, according to Eq. (|245p : 



TpiUi <E) U2) 



(279) 



k 

E 



^C0S2 % 



,2 9k 



n2 Xfc 



\^sin^ ^ cos2 V^in ¥ 

'cos2fpi?+sin2fp(f| /cos 
3in^^p<^;)+cosfp(^2V V^in^ 



„2 Xfc 



2 Xfc TrC^) 

"2"^'21 

2 Xfe 

P21 



sin^ 2aLp<2^ 

- cos ^p^2? , 
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which eventually leads to the following equation 



1/ o di 2 ^2 2^1 9 ^2 /I /I / \\ 

- I cos — sm h sni — cos h smfc'i sin6^2 cos[tpi — (p2) ] 

E\ ( 2^k~{k) , . 2^k~{k)\/ 2Xk~(k) , . 2 Xfe~<fc)\ .„cn^ 
Afe (^cos — pi^^+sm YPi2 j (cos — p^2^^ + sm ^^22 ) (280) 

with an infinite number of unknown variables Xk,dk,Xk,Pii , P12 j P21 ' P22 ■ This 
equation has no solutions, so the decomposition of Eq. (|279p is impossible for 
this tomographic probability vector: the state is entangled. The proof is given 
in the next subsection. 

To resume, we have presented examples of tomograms for two pure simply 
separable states. One state corresponds to both spins directed along z-axes. The 
density matrix of this state in the natural basis has the form of Eq. (|247p and 
the tomogram is given in the form of probability vector (j253[) with the tensor 
product form of probability vectors of Eq. (|254p . Another state corresponds to 
both spins directed along the x-axis. The density matrix of this state has the 
form (|256p and the tomogram is given by (|263p . The tensor product form of 
the tomogram is given by Ea. (|264p . The separable, but not simply separable, 
mixed state with density matrix (|265p has the tomogram given by Eq. (|266p in 
the form of convex series (j245p . with only two nonzero terms in the series. 

An example of entangled state is given by the pure state (|267p . with density 
matrix (|268p and tomogram (|279p . 

4.6 Inequalities 

To prove that equation (|280p has no solution let us use the following inequality 
valid for stochastic 4 x 4-matrices which are a tensor product of two stochastic 
2 X 2-matrices, i.e. 

(281) 

Let us introduce the matrix 



M = 


Ml € 


iM2. 




/I 


-1 


-1 




1 


-1 


-1 


1 


1 


-1 


-1 


1 


V-1 


1 


1 


-V 



Then we prove that 



(282) 



|Tr(M/o)| < 2 (283) 



is a necessary condition for the validity of Eq. (|28ip . 
In fact, the stochastic 4 x 4-matrix has the form: 



P1M2 qiA'h 



M = '''Z , (284) 

U2M2 (72M2 ' ' ^ ' 
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where the stochastic matrices Mi and M2 read 



M,^(P' '');M,^h lA. (285) 

\P2 92/ \S2 t2j 

The columns of Mi and M2 are probabihty vectors: pi +p2 = qi+q2 = S1 + S2 = 
ti+t2 = 1 and aU the matrix elements are non- negative. The trace in Eq. (|283p 
reads: 

Tr(Af/o) = Piisi ^ S2) - P2{si - S2) + piih - t2) - P2{h - t2) (286) 
+ qi{si - S2) - 92(51 - S2) - qiih - <2) + q2{ti - <2) 
= (Pi -P2)[(si - S2) + (ti - t2)] + {qi - 92) [(si - S2) - (tl " i2)] 
=: p{s + t) + q{s - t) 

The differences p, q, s, t of the probabilistic distributions satisfy respectively the 
inequalities 

\pi -P2I < 1 , ki - 92I < 1 , |si - S2I < 1 , \ti -t2\<l (287) 

To prove that the modulus of the sum in Eq. (|286p does not exceed the 
number 2, consider the function /: 

f{p,q,s,t):^p{s + t)+q{s~t) (288) 

that is a harmonic function of the four variables p, q, s, i, which arc constrained 
to belong to the hypercube K/^ — {\p\ < 1, \q\ < 1, \s\ < 1, \t\ < 1} . Then max- 
imum and minimum values of / on are reached on the boundary of the 
hypercube. Note that / has the special properties that, when any number of 
its variables is taken constant, / is still harmonic in the remaining variables. 
Therefore / is harmonic when restricted on each face of the boundary of the 
hypercube, which is an hypercube K3 of one less dimension. In each face K3 the 
max and min of / lie on its boundary. By repeating this argument, eventually 
the max and min of / are found to lie on the vertices of the initial hypercube 
K4, which are 2^ = 16 points. For the given function it is then a trivial matter 
to check that max / — 2 and min / = — 2. 

The argument has an immediate generalization to any number of dimensions: 
Proposition: Let / be a function of n variables (a;i, ...,a:„) harmonic on the 
hypercube Kn — {\xi\ < l;i = l,..,rt}, such that / is still harmonic when 
restricted to any face A's,0 < s < n — 1, belonging to the boundary of Kn. 
Then / has its extrema on the 2" vertices of Kn- 

^From inequality (|283p it follows that 



Tr 



(/o XkMk) I < 2 (289) 



where Xk > OjJ^k^^k — 1 and the matrices have the form of Eq. (|28ip . 
Then, the proof of non-existence of solution to Eq. (|280p is reduced to check the 
violation of inequality Eq. ()289p . 
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Let us construct the 4 x 4- matrix M using the probabihty vector (|279p . We 
take this vector for four pairs of argument matrices, namely: 



Thus, the constructed matrix M is a function of eight angles: the angles 0a, <^a 
which are the Euler angles determining the 2 x 2- matrix u'\^\ those 9d,^d 
determining the 2 x 2-matrix u['^\ while Ob, (fib determine the 2 x 2-matrix 1/2^^ 

('2) 

and 6'c, ^Pc the 2 x 2-matrix C/j • 

If one assumes that the equality (|279p is valid, then the constructed matrix 
(|290p has the form of a convex sum J2k ^kMk satisfying the inequality (|289p for 
all values of the eight angles 9a,'Pa,0d,'Pd,db,(pb,Oc,ipc- On the other hand we 
know the explicit form of such a matrix. In fact, the first column of this matrix 
is 

f Xab \ 



r(2) 



(2) 



^(2)^ 



V 



where 



Xab - 2 



2 f 6 



(291) 



3bCos{^a-Vb) , (292) 



the second column is obtained from the first one (|29ip by the following replace- 
ment: 

a , b^c , (293) 
a-^ d , b^b , (294) 



the third column by: 
and the fourth column by: 



(295) 



Taking the trace of the obtained matrix M with the matrix Iq of (|282p . we 
get the following function of the eight angles: 



B = A{Xab + Xac + Xdb ~ Xdc) - 2 



(296) 



Now we look for values of the angles 9a, ^a, 9d,^d, 9b, ^b, 9c, for which \B\ 
exceeds 2 (in fact |82] the maximum of this function is equal to 2\/2). Such 
values do exist and this implies that the hypothesis ()279|) for the given prob- 
ability vector is false. In particular, the maximum of \B\ is achieved when 
Xab = Xac = Xdb = 1 ~ Xdc, and the Corresponding stochastic matrix reads 



M = 



2+^2 


2+^2 


2- 


fy2 


2-^2 


8 

2-V2 


8 

2-\/2 


2- 


8 

-V2 


8 

2+V2 


8 

2-\/2 


8 

2-\/2 


2^ 


8 

-V2 


8 

2+V2 


8 

2+^2 


8 

2+\/2 


2- 


8 

f\/2 


8 

2-\/2 



(297) 
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It means that Eq. (|280p has no solution. The incquahty 

\B\ < 2 (298) 
is called Bell inequality [53] (or CHSH inequaUty 84J). 
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5 Conclusions 



To conclude we point out the main aspects discussed in the paper. We reviewed 
the probability representation of quantum states, in which wave functions or 
density states are replaced by tomographic probability distributions containing a 
complete information on quantum states. The mathematical mechanism to con- 
struct all the possible probability descriptions of quantum states was clarified. 
It amounts to constructing complete sets of rank-one projectors in the Hilbert 
space of operators acting on the underlying Hilbert space of state vectors. These 
sets are complete or overcomplete. The tomograms depending on continuous 
variables were shown to be used both in the classical and quantum domains. 
The set of all tomographic functions describing the quantum states, the classical 
states or no states at all was characterized. The characterization is expressed 
in terms of inequalities, given by Eq.s ([73|), l|74)) . We suggested a method of 
direct experimental checking the Heisenberg uncertainty relations [7T]. In ex- 
periments [HI [T51 [751 the optical tomogram of a photon state T{X, 9) was 
measured in order to find the Wigner function yV(q,p) by means of the Radon 
anti-transform of the tomogram. We pointed out that the new inequality (|157p 
for the dispersion of the homodyne photon quadrature X — cos 9Q + sin 6'P, ex- 
pressed in terms of integrals containing the directly measured optical tomogram 
T{X, 9), provides the method of checking Heisenberg uncertainty relations. The 
accuracy of the experiment on checking the uncertainty relations is in one-to- 
one correspondence with measuring optical tomogram in the above experiments. 
The inequality (|157p can be violated if one uses tomograms of classical states. 
The Heisenberg uncertainty relation checking can be possibly accompanied by 
a checking of the inequality for tomographic entropy discussed in ! 85j . 

The entanglement of spin system states was given in terms of properties of 
the spin tomograms. The Bell inequalities were shown to reflect the properties 
of entanglement in terms of properties of joint probability distributions (spin 
tomograms) to be expressed (or not expressed) as convex sums of joint proba- 
bility distributions without correlations. The quantum spin states provide some 
bounds for the correlations. The Cirelson [82] bound 2\/2 corresponds to the 
Bell number characterizing a maximally entangled state of two qubits expressed 
in terms of the system tomograms properties. The tomographic-like joint prob- 
ability distributions for which the Bell number is greater than that bound (the 
maximum can be equal to 4, see [86]) do not correspond neither to states with 
quantum spin correlations nor to classical correlations, characterized by a bound 
equal to 2. The stochastic matrix which has the Bell number 4 is, for example, 
of the form: 

i 

I 

Vo ! 1 oy 

Suppose, in an experiment on measuring the spin state tomogram of two qubits, 
to get the probability vector (|25ip and the stochastic matrix, whose columns are 
the values of this vector for 4 pairs of Euler angles, to be just the above matrix 



M 



(299) 
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M. This result contradicts quantum mechanics. Moreover, the 'classical' states 
of two qubits have the bound 2. Thus, the value 4 of the Bell number corresponds 
to a tomogram which is neither quantum nor classical, and this is similar to the 
example with the scaled Wigner function for continuous variables discussed in 
section [23] We mean 'classical' in the following sense: in classical probability 
theory, given a distribution function of two random classical variables, there are 
two possibilities. One is that the random variables are uncorrelated, i.e, the joint 
probability distribution has a factorized form as a product of two probability 
distributions, describing statistical properties of each of the random variables. 
Another possibility is that the random variables are correlated, so the joint 
probability distribution is a convex sum of factorized probability distributions. 
The correlations described by such sums we call 'classical'. In a sense, it is 
a terminology that refers to a class of situations in which the joint probability 
distribution of two correlated random variables is thought of as a convex mixture 
of distributione without correlations. 

This picture for discrete spin variables is analogous to the picture with con- 
tinuous variables, where tomograms violating both inequalities ([75)) . ([7^ do not 
correspond neither to classical nor to quantum states. 

We have shown that in tomographic approach to quantum mechanics an im- 
portant role is played by semigroups and their orbits in simplexes, since from 
a geometrical point of view the tomograms, being probability distributions, are 
points of simplexes and their dependence on extra parameters provides some 
domains in the simplexes as orbits of semigroups. For qubit systems the semi- 
groups are obtained from unitary matrices taking the map of their elements 
onto their square moduli. There exist generalizations of symplectic tomographic 
maps both for classical mechanics of a top |87j and quantum mechanics with 
using maps with curvilinear coordinate lines Radon transforms |H5] . The analy- 
sis presented in this paper may be extended to those cases of generalized Radon 
transform. We hope to study these maps and their properties in future papers. 

Among the variety of open problems, we should mention few of them which 
are more crucial for the full equivalence of the tomographic picture of quantum 
mechanics with the existing ones. We do not have yet a complete and au- 
tonomous characterization of tomograms whose inverse Radon transform gives 
rise to a Wigner function or to a classical probability distribution. Similarly 
we do not have a complete characterization of the continuity property of the 
Radon transform and its inverse. Some of those aspects play a relevant role to 
establish the topological properties of the star-product among the observable 
functions and the action of the observables on the states. We shall consider 
these interesting problems in a forthcoming paper. 
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